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Abstract 

We study the dynamics of a polymer or a D-dimensional elastic manifold dif- 
fusing and convected in a non-potential static random flow (the "randomly driven 
polymer model"). We find that short-range (SR) disorder is relevant for d < 4 for 
directed polymers (each monomer sees a different flow) and for d < 6 for isotropic 
polymers (each monomer sees the same flow) and more generally for d < dc{D) in 
the case of a manifold. This leads to new large scale behavior, which we analyze 
using field theoretical methods. We show that all divergences can be absorbed in 
multifocal counter-terms which we compute to one loop order. We obtain the non 
trivial roughness C,, dynamical z and transport exponents in a dimensional ex- 
pansion. For directed polymers we find C, ~ 0.63 (d = 3), C ~ 0.8 [d = 2) and 
for isotropic polymers C, ~ 0.8 (d = 3). In all cases z > 2 and the velocity versus 
applied force characteristics is sublinear, i.e. at small forces v{f) ~ with cf) > 1. 
It indicates that this new state is glassy, with dynamically generated barriers lead- 
ing to trapping, even by a divergenceless (transversal) flow. For random flows with 
long-range (LR) correlations, we find continuously varying exponents with the ratio 
dh/dT of potential to transversal disorder, and interesting crossover phenomena be- 
tween LR and SR behavior. For isotropic polymers new effects (e.g. a sign change of 
C ~ Co) result from the competition between localization and stretching by the flow. 
In contrast to purely potential disorder, where the dynamics gets frozen, here the 
dynamical exponent z is not much larger than 2, making it easily accessible by simu- 
lations. The phenomenon of pinning by transversal disorder is further demonstrated 
using a two monomer "dumbbell" toy model. 
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Figure 1.1: elastic manifolds (polymers £> = 1) in random flows: (a) directed polymer (b) isotropic chain. 



In this paper, we study the model of a polymer diffusing in a static random flow and its 
extension to an elastic manifold of arbitrary internal dimension D. The velocity pattern 
of the flow v{r) is constant in time and leads to convection of the polymer in addition 
to diffusion. We are interested in the general case of a non-potential flow: the extreme 
example is the hydrodynamic divergenceless flow, with V • v{r) — 0, but mixtures of 
potential and divergenceless flows are also considered. See figures 1.2, 1.3 and 1.4 for a 
visualization. The physical motivation to study this model, which we call a "randomly 
driven polymer", are the following: 

First, it is a prototype of an out of equilibrium system in the presence of quenched 
disorder which generalizes the much studied problem of an elastic manifold in a random 
potential [1-11], i.e. for the present model the case of a purely potential Row . This problem 
has been very popular in the last decade due to its numerous physical realizations such as 
vortex lines in superconductors, interfaces in random magnets and charge density waves 
[12-15]. Equally fascinating, albeit more formal, relations exist to growth processes and 
potential fluid dynamics, via the Cole-Hopf transformation onto the Burgers and Kardar- 
Parisi-Zhang equations [16]. In addition, it is one of the simplest models which exhibits 
glassy properties, typically ultra slow dynamics, anomalously small response to external 
perturbations and pinning [2,12,17], which are also present in more complex systems, such 
as spin glasses [18-20] or random field systems [21,22]. However little is known about what 
of the glassy properties remains in the presence of non-potential disorder. This question 
is important for numerous physical systems such as driven systems in the presence of 
quenched disorder [23-26] or domain growth in the presence of shear [27]. The model 
studied in this paper is simple enough to allow for some answers to this question. 

The second motivation concerns real systems of elastic objects convected by hydro- 
dynamic flows. Recently, experiments on polymers in spatially varying stationary flows. 



1 Introduction 
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Figure 1.2: Longitudinal (potential) disorder, obtained from the disorder of Figure 1.4 through application 
of the longitudinal projector. 

such as elongational flows have been performed [28-31] and these experiments can be 
extended to longer objects and more complex flows. Of course the model studied here 
is a toy model in the sense that it disregards hydrodynamic forces (the back reaction of 
the polymer on the flow) and additional time dependent components of the flow. The 
main interest is to illustrate in a simple way some of the fascinating properties arising 
from the competition between the convccting flow and the elasticity of the structure. We 
believe that some of these properties, for instance the existence of preferred regions in the 
flow, will be present in more realistic situations. Experimentally, inhomogeneities in the 
polymer distribution and non-linearities in the response of the drift velocity v{f) to an 
uniform applied external force / arc worth being investigated [32]. 

In the case of nonpotcntial disorder new effects arise from the competition between 
on one hand disorder and elasticity - which tend to create pinned or frozen states - and 
on the other hand the energy which is constantly pumped into the system and which 
tends to destroy the glassy properties. Very few systematic studies of these effects exist 
at present. One extensive work was performed in the context of driven lattices [25.33]: it 
was shown that some glassy features, such as pinning in the displacement transverse to the 
motion, barriers and strongly non linear transverse velocity force-characteristics, survive 
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Figure 1.3: Transversal disorder, obtained from the disorder of Figure 1.4 through appUcation of the 
transversal projector. 

in the moving state. Some glassy features were also found to survive in infinite range 
random spin systems with non relaxational dynamics [34]. Another example, studied 
long ago [35-43] and recently solved analytically in the limit of d — > oo for long-range 
(LR) disorder [44], is the problem of a single particle diffusing in a quenched random flow, 
the D = limit of the present model. This problem was found to exhibit some glassy 
features, also seen in recent simulations of a related spin model [34]. In the present paper 
we generalize the single particle study to an elastic manifold, which is expected to lead 
to new physics due to the crucial role of internal elasticity in generating barriers. It also 
generahzes the studies of the dynamics of self-avoiding manifolds [45,46] to a quenched 
disorder situation, though self-avoidance is found to be perturbatively less relevant in 
most of the regimes studied here. 

In this paper we work directly in finite dimension using dynamical renormalization 
group techniques, and establish the existence of new universal fixed points for short- 
range disorder. We also treat the case of long-range disorder, which leads to a fine of 
fixed points. In all cases we compute explicitly the non trivial roughness exponent, the 
dynamical exponent and the velocity force characteristics at small applied force /. We 
study two cases: either the polymer is directed (each monomer sees a different fiow) 



5 




Figure 1.4: Isotropic disorder, obtained through superposition of 7 x 7 Fourier-modes. One observes 
purely transversal ("turbulent") regions and potential wells. 

or it is isotropic (all monomers are in the same flow) as depicted on Fig. 1.1. A short 
account of our results has already appeared in [48], which was presented using a Wilson 
RG technique as in related studies [33,49]. The aim of the present paper is to give a 
detailed account of our results using a field theoretical derivation. This derivation, based 
on the powerful multilocal operator product expansion (MOPE) techniques developed in 
the context of self- avoiding manifolds [50-54,45,46], deserves a detailed presentation, as 
it is generally useful to treat problems involving manifolds. In addition, we give here a 
detailed treatment of the crossover between long-range and short-range disorder, as well 
as an in depth discussion of the results. 

This article is organized as follows: In section 2 we define the model, the physical 
observables and present a simple analysis based on dimensional and Flory arguments. 
Section 3 is devoted to the field theoretic treatment of the model. We start by defining 
the generating functional (section 3.1), and then show that divergences are short-range 
(section 3.2). This allows to identify local (section 3.3) and bilocal (section 3.4) counter- 
terms, summarized in section 3.6. Their divergent parts are calculated in sections 3.7 
and 3.8, followed by a discussion of the particle limit in section 3.9. This yields the 
RG-equations in section 3.10 and 3.11. 
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The section which contains a presentation and a discussion of the results, and can be 
read independently is section 4. There we analyze the RG equations, report the results 
of our calculation and specify them to the four cases of interest (directed and isotropic 
disorder, short-range and long-range correlations). The crossover from long-range to short- 
range correlated disorder is discussed in section 4.6. We compare with previous results 
for Hartree limit [d = oo) and particle case {D = 0) and discuss the crossover from 
short-range to long-range correlated disorder. This is followed by an analysis of the 
effect of self-avoidance in section 4.7. The section terminates with the study of a toy 
model, involving 2 particles in a (directed) random flow which illustrates how barriers are 
generated dynamically. Conclusions and perspectives are presented in section 5. Some 
technicalities are relegated to appendices A, B and C. 

2 Model and physical observables 
2.1 The model 

We consider a manifold of internal dimension D parameterized by a d-component field 
ra{x). The polymer corresponds to D = 1 {x labels the monomers), and a single particle 
to D — 0. Membranes have D — 2. We study the Langevin dynamics (for the dimensional 
analysis see section 2.3) 

Taix, t)^X (Ar^ix, t) + rjaix, t) + F^[r{x, t),x]) , (2.1) 

where the term proportional to Ara{x,t) is the derivative of the internal (entropic) elas- 
ticity ("stiffness"). The Gaussian thermal noise 7]a{x,t) satisfies 

{rj^ix, t)r)f,{x', t')) = j5^f,5{t - t')S''{x - x') . (2.2) 

Angular brackets denote thermal averages and overbars disorder averages. 
FQ,[r, x] is a Gaussian quenched random force field of correlations: 

F^[r,x]Fp[r',x'] = A^^ir - r')h{x - x') . (2.3) 

We study two main cases of interest illustrated in Fig. 1.1. If the manifold is directed 
(e.g. a polymer oriented by an external field), then 

h{x - x') = S^{x - x') . (2.4) 

If the manifold is isotropic, (e.g. a Gaussian chain in a static flow), the force fleld does 
not depend on the internal coordinate 

F^[r,x]=Fa[r] (2.5) 

and 

h{x-x') = l. (2.6) 

Note that the isotropic polymer (or more generally the manifold) in a static flow can also 
be seen as a directed polymer in a flow completely correlated (i.e. constant) along its 
internal dimension. 
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We consider a statistically rotationally invariant force field with both a potential (L) 
( "electric" ) and a divergence- free (T) part ( "magnetic" ) which depends only on the dis- 
tance \r — r'\. Both parts contribute separately to the correlator: 

A.Mr) = -drAA^{\r\) - {5a,Y.dl^ ' drA,)^^m ■ (2.7) 

7 

In Fourier space we use 

A^p{r) = / A«^(A;)e^^"- (2.8) 

Jk 

where we specify the normalization of the fc-integral below. There are several cases of 
interest for the correlation of the random force. 

Model SRF: force with short-range correlations 

The force correlator scales like a 5-distribution Ao/3(r — r') ~ 5'^{r — r'), with however 
a non-trivial index structure. In Fourier space the correlator reads at small k 

^ap{k) = g^Plp{k) + gT:Plp{k) , (2.9) 

where Papik) = kakjs/k'^ is the longitudinal and -P^(A;) = 5^/3 — fapi^) transverse 
projector. As we shall show below, this is the generic situation, since it is relevant for all 
short-ranged correlations, even those, for which one starts with a force correlator which 
is formally shorter ranged than the 5-distribution (e.g. decaying faster than l/r*^). 

Model SRP: potentials with short-range correlations 

In this case, the |Jotentza/ correlator scales like a ^-distribution Aj^^T{r—r') ~ S'^{r — r'). 
The corresponding force-correlator reads in Fourier space at small k 

Aap{k) = gLkakp + MSapk'' - kak^) = e [gLPapik) + gTPapik)] . (2.10) 

Except for = 0, which is preserved by a fluctuation dissipation theorem [12,5,55], 
this model is unstable and will renormalize to model SRF discussed above (except in the 
transversal case gh = for the particle D = 0). 

Model LR: forces with long range correlations 

In this model, correlations at large scale obey a power law Aa/sir — r') ~ |r — r'|~", 
with a < d. The Fourier transforms of the correlators with respect to r behave, for small 
values of their argument k, as 

A^^ik) ~ \kr'{g^Pt^{k) + grPlpik)) . (2.11) 

Let us also note that in order to simplify calculations, we normalize the integral over 
k such that 



Ik 



2.2 Physical observables 

We study the time translational invariant steady state. It can be constructed by taking 
a finite size sample L and considering the large time limit first, before taking the system 
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size to oo. We will not address here the subtle issue of whether, or in which case, these 
limits commute (for a partial discussion in the case > see e.g. [44]). The radius of 
gyration (or roughness) exponent C is then defined by 



{{r{x, t) - r{x', i))2) ~ |x - ^ , (2.13) 
and the single monomer diffusion exponent u by 



{{r{x, t) - r{x, f ))2) r.\t- t'l^" . (2.14) 

We assume scahng behavior 



{{r{x, t) - r(0, 0))2) ^ B{t/x') (2.15) 

with ( — zu. 

We also consider the anomalous dimension of the elasticity of the membrane, /3 (for 
the precise definition sec Eq. (3.108)), which evaluates to in the directed case. 
The drift velocity under a small additional applied force / in Eq. (2.1) is 



V {r{x,t))/t . (2.16) 

We find 

vr^f (2.17) 

at small /, with 

>1, (2.18) 



^ 2-C + /3 
indicating trapping of the polymer by the flow. 

2.3 Dimensional analysis 

Power counting at the Gaussian fixed point without disorder yields (denoting e.g. [r]^ the 
dimension of r in units of x) 



[rl = Co = , [r], = 1^0 = , Co = = [rX = — , /? = 

(2.19) 

for D < 2. The dimensions of the couplings are 

£^^^(D, d) = - = - [^T]f = 2 - Coa (2.20) 

6'-(D, d) = - = - [^T]f = 2 + - Coa (2.21) 

in the directed and isotropic case respectively. Thus short range disorder (formally a :— d) 
is relevant for positive £, i.e. when (see Fig. 2.1) d < d'^^^{D) for directed manifolds and 
d < (4f°(-D) for isotropic manifolds, with: 

4 4 _|_ 2D 

dTiD) = ^— ^ , d-^{D) = . (2.22) 
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Figure 2.1: Regions in the {d,D) plane where SR disorder is relevant for a directed manifold. Inset: 
Regions in the (a, d) plane where SR or LR disorder is relevant for the directed polymer D = 1. 



Thus we find that short-range disorder is relevant for directed polymers when d < 4 and 
for isotropic polymers for d < Q. This is interesting because it means that in physical 
dimensions {d = 2, d = 3) polymers should be affected by quenched random flows while 
for particles short-range disorder is irrelevant for d = 3 and only marginal in c? = 2. 
Disorder can also be relevant if the interaction is long enough ranged: disorder is relevant 
for a < ac{D) — dc{D) for LR disorder. 

A qualitative idea of the effect of disorder is obtained by using Flory-like estimates. 
Balancing in Eq. (2.1) friction or elasticity with disorder ^ yields 

(dir ^ 2vf = ^ , Cf° = 24'° = tA-. (2-23) 

for directed and isotropic manifolds, respectively. Results for LR disorder are obtained 
upon replacing d ^ a. Note that this argument is too rough to distinguish between the 
potential or divergenceless nature of the flow. It suggests however that the manifold will 
be stretched by disorder. As we will see below, the precise calculation will indeed confirm 
that the manifolds are stretched in most (but not all!) of the cases, and will yield the 
deviations away from Flory. 

The expected "phase diagram" in the d, D plane is illustrated in Fig. 2.1 for directed 
manifolds. The new glassy regime extends from the line d = dc{D) to the upper critical 
(internal) dimension D = = 4 beyond which disorder is perturbatively irrelevant. 

^The idea of the Flory-estimate is to suppose that all terms in the equation of motion scale in the same 
way, and by this to obtain the critical exponents. For directed manifolds in short-range correlated forces, 
one sets dtr ~ Ar ~ F{r,x), replaces dt by 1/t, A by l/x"^ and from Eq. (2.3) F{r,x) by r~'^/'^x~^/'^ to 
obtain r/t ^ rjx^ ~ r~'^/'^x~^/'^ , and the exponents stated in the main text. For isotropic manifolds in 
short-raugo correlated forces, the same argument reads r/t<^r/x'^'^ r~'^l'^ . 
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For isotropic manifolds the diagram has the same topology except that the upper critical 
internal dimension is pushed to = +00. We know of three ways to investigate per- 
turbatively the nature of the glassy regime in Fig. 2.1. First one can start from large d 
(any D), which was done in [44]. The second approach is to start from D — A using a 
functional RG method [73]. The last approach is to start from the line d — dc{D) which 
is the aim of this paper. 

3 Field theoretic treatment of the renormalization group equations 

3.1 Definition of the field-tlieoretic model 

We start from the equation of motion (2.1) and convert it into an effective field theory 
following [56]. The effective dynamical action in prepoint (Ito) discretization reads 

J[r,f]^ [ faix,t){r^{x,t) + X{-Araix,t)+r]aix,t) + F^[r{x,t),x])) . (3.1) 

The noise r]{x^t) as well as the disorder force F{r) are Gaussian, see Eqs. (2.2) and (2.3), 
and can thus be integrated over. This gives the dynamic functional 

J[r,f]= / f{x,t){f{x,t) -X/^r{x,t)) - \f{x,tf 

Jx,t 

-\\^ j f«(x, t) A,^(r(x, t) - r{y, t'))h{x - y)fp{y, t') . (3.2) 

^ Jx,y,t,t' 

Here and in the following, contraction over indices is implied, where ever confusion is 
impossible. In order to simplify notations, we will in the following suppress the factor of 
A, i.e. set 

Xt^t . (3.3) 

This is not problematic, as A always appears with time. At the end of the calculations 
one has to replace t by At which is necessary in order to get the renormalization factors 
correct. Note also that the factor | disappears, if we use the time-ordered vertex. 
The disorder correlator is 

A«Mr(x, t) - r(y, t)) = e^Kri^'^^-rM) ^g^P^^(k) + g^Plpik)) \kr''h{x - y) . (3.4) 

The function h{x) may be specified to the directed case with 

h{x - x') = 5^{x - x') (3.5) 

or to the isotropic case 

h{x - x') = 1 . (3.6) 
Eventually, one may also want to study an interpolating case, where 

h(x - x') = 5(x± - x'^) . (3.7) 
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In order to simplify notations, we introduce the following graphical symbols 

— f(x, t)r{x, t) 
"'^'^ — f{x, tY 

fa(x,t)e^'=('^(=^'*)-"(^'*'))f^(y,t')|A;r"'^/i(a: - y) (3.8) 

Jk 

P^p{k)fa{x, t)e^'=('-(^'*)-'-(^'*'«f^(|/, t')\kr'h{x - y) 
[ P,T^(A;)f,,(x,i)e^'=(^(-'*)-'-(^'*')V>(y,i')|^r-'/i(^-y) ■ 

Jk 

The action takes the symbolic form 

Jx,t 2 Jx,y,t,t' ^ ' 

The free theory ((^l = S't = 0) is easily investigated to obtain the free response R{x^ t) 
and correlation- function C{x,t): 

R{x,t) = ^(r(a;,t)f(0,0))„ (3.10) 
C{x,t)^^(^\{r{x,t)-r{0,0)r)^ , (3.11) 

where (. . .)q denotes averages in the free theory. Explicit analytical expressions are given 
in appendix A. Here it is only important to note that for < D < 2 both functions 
are well defined in the continuum. C{x,t) vanishes ior x — t — 0, and is monotonically 
increasing in both arguments. They are non-trivially related by the fluctuation-dissipation 
theorem 

e{t)^^C{x,t) ^ R{x,t) . (3.12) 

3.2 General considerations about renormalization 

Before embarking upon calculations, let us follow [51,45] to show that divergences are 
always short-range. This is necessary for the renormalization program to work. To this 
aim write down the perturbation-expansion, the interaction vertex being 

:= / r{x, t) eik{r{x,t)-r{y,t')) ~^y^ ^ (g^^^g) 
J k 
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where for the sake of transparency, we have neglected the index structure. The reader 
will easily be able to generalize the analysis given below to the model defined in Eq. (3.2). 
The perturbative expansion of an observable O can then be written as 

(e)) = Norm^^/ (O— --")o , (3-14) 
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where (. . .) denotes the expectation vahie with respect to the free, non-interacting theory, 
and the normahzation Norm has to be chosen so that (1) = 1. The integral is taken over 
all arguments of the interaction vertex. We claim that divergences only occur at short 
distances and short times. To prove this look at a typical expectation value with fixed Xi 
and ti 

{O^ ^«)^ = ^ / f^^xi-x^^ti-t^^k^kje-^^.^"^'^''"^ , (3.15) 

where each contribution consists of a function /„, which is a polynomial product of cor- 
relation and response functions and fc's and an exponential factor, with 

Qij = C{xi Xj,ti tj) . (3.16) 

As fa is a regular function of the distances, divergences in the /c-integral in Eq. (3.15) can 
only occur if Qij is not a positive form. We will show that Qij is a positive form for all ki 
which satisfy the constraint 

E^i = 0- (3.17) 

i 

This constraint always holds, see Eq. (3.13). For equal times it is just the statement that 
the Coulomb energy of a globally neutral assembly of charges is positive by identifying C 
with the Coulomb-propagator and ki with the charges. In the dynamic case, write 

Q = f / ^ 1 f^p{x,-X,)+iu.(t.-t,) _ ^\ (3_;L8) 

J (27r)^ J 27r a;2 + (p2)2 V J ^ ^ 

The exponential in Eq. (3.15) then reads 

Vfc,fc,g,, = / / — — '^—^y^k.k, (e^p(-.-.)+-(*^-%) _ i) 

r d^p r dcu 2 
^ J (27r)^ J 2^^7TW 

To get the second line, Eq. (3.17) has been used. Note that again due to Eq. (3.17), the 
integral is ultraviolet convergent and thus positive. It vanishes if and only if the charge- 
density, regarded as a function of space and time, vanishes. This is possible if and only if 
endpoints of the dipoles (which form the interaction) are at the same point in space and 
time. 

Thus when integrating in Eq. (3.14) over Xi and ti, no divergence can occur at finite 
distances. To renormahze the theory, only short distance divergences have to be removed 
by adding appropriate counter-terms. Thay are analyzed via a multilocal operator product 
expansion (MOPE) [50,51]. For a detailed discussion of the MOPE and examples see 
[50-52,54,47]. The main result of this analysis is, like in scalar field theory, that in 
any order of perturbation theory, only a finite number of counter-terms is necessary and 
that these are exactly those terms, which are already present in the original theory. The 
theory is thus multiplicatively renormalizable. To actually evaluate the counter-terms, a 



(3.19) 
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multilocal operator product expansion is used, which is described in the next sections. The 
general idea is always to study the behavior of n disorder-vertices on approaching their 
end-points. The resulting term is expanded into a product of a (multi-local) operator, 
times a MOPE-coefficient, which is a homogenous function of the distances involved. 
Power-counting then indicates that as long as the operator is relevant or marginal in the 
RG-scnsc, the integral over the distances possesses an UV-divergence, and thus requires 
a counter-term in the Hamiltonian. On the other hand, when the operator is irrelevant 
in the RG-sense, then the integral over the MOPE-coefficient is globally UV-convergent, 
and no counter-term is needed. This procedure is constructed such, that possible sub- 
divergences are taken care of recursively as in the BPHZ-scheme. A formal proof of this 
statement is given in the context of a simpler model in [51], see also [47]. 

3.3 Divergences associated with local operators 

We now analyze our specific model, using the techniques of the multilocal operator product 
expansion (MOPE) [50,51]. The first class of divergences stems from configurations where 
the two end-points of the interaction vertex are approached. The interaction vertex is: 

..^^L^ f f,(x,t)e^'=['-(^'*)-^(^'*')lf^(|/,t') X 

Jk 

X {Pap{k)9L + Plp{k)gT) Ikr'hix - y) . (3.20) 

For short-range forces set a — d. 

In order to extract the divergences for small x — y and t — t', the first possibility is not 
to contract any response-field. We then start by normal-ordering the r.h.s. of Eq. (3.20). 
As usual in statistical field theory, we define the normal-ordered operator : O : starting 
from the operator O by 

:0: — O — all self-contractions (3-21) 

where the self-contractions are taken with respect to the free theory. The idea is that in 
order to calculate the insertion of :0 : into some expectation value only contractions of 
: O : with the other operators in this expectation value have to be taken. This implies 
that the operator : Oi{x, t)O2{0, 0) : is always free of divergences for {x, t) —>■ (0, 0). 
Using the identity 

^ikr{x,t) ^-ikr{y,t') _ . ^ikr(x,t)^-ikr{y,t') . ^-k'^C(x-y,t-t') (3 22) 

one can expand the normal-ordered vertex-operators on the r.h.s. for small x — y and t — t'. 
The leading contribution is 

le-k'c(x-y,t~t') ^ (3 23) 

yielding the first term in the short-distance expansion of Eq. (3.20) (for the normalization 
of the A;-integral cf. Eq. (2.12)): 

/ e-^^^(-^'*"*') ■M^,t)ffi{y:t'): {PUk)9L + Plf3{k)gT)\krM^-y) 

Jk 



k 

Jx + y t + t'^^ 
= r — - — , — 



9t 



^(/t ~ ^) ~'~ ^^^)^^^ ^ ^ t') "'^'^h{x — y) + subleading terms 
1 - + ^L^) C{x -y,t- t'y^'^h{x -y) + subleading terms . (3.24) 
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By going from the first to the second hne, two things have been used: First, since the 
integral over k is rotationally invariant, only terms with a = (3 survive, with a geometric 
factor of 1/d for the transversal and (1 — 1/d) for the longitudinal projector. Second, 
since :f(x, t)f{y, t') : is normal- ordered, i.e. all contractions between the two fields f have 
been eliminated and no short-distance singularities remain, it can be written as 

w .n Jx + y t + t'V Jx + y t + t'\ . .Jx + y t + t'\ 

■■r{x,t)r{y,t): = :r ——,—— :+ :r [2(x - y)V] r : 



^(x + y t + t'" 



+ 0{{x-yY)+0{{t-t'Y)+0{x-y)0{t-t') (3.25) 



From this expansion, only the first term has to be retained and will give a correction to the 
marginal operator = in the action. The following terms are irrelvant (in the RG- 
sense) with respect to this dominant term and can therefore be neglected. Dropping the 
explicit normal ordering of which is a pure number, and zero by analytic continuation, 
leads to Eq. (3.24). 

Note that when expanding the product of two operators as sum of normal ordered 
ones, this sum runs over all possible contractions. The second contribution to the normal- 
ordered product of Eq. (3.20) is therefore obtained upon contracting one response field. 
Due to causality, this must be the field with the smaller time-argument. For simplicity, 
let us take t > and put y — t' — 0. By the same procedure as above, we obtain the 
following contribution: 

/ :f«(a;,t)e^'=(^(^'*)-'-('^'°»: i?(a;, t)(zfc)^e-^'^(-'*) (P,L^(fc)5L + Pj^(fc)5T)|fcr"'^/^(x - y) . 

(3.26) 

The next step is to expand :e^KH^,t)-Ho,o)) ■ about (x, t). (It is important to expand about 
{x, t) as otherwise f(x, t) had to be expanded, too.) This expansion is 

:eik{r(x,t)-rm).^l + ^if^^_^ (^tf^{x,t) + {xV)r^{x,t)-^{xV)\{x,t)^+suhlesidmg terms . 

(3.27) 

Upon inserting Eq. (3.27) into Eq. (3.26) and integration over k, only terms even in k 
survive. We can also neglect the term linear in x, which is odd under space reflection. 
The remaining terms are 

^ ■.fa{x,t) (^tr^{x,t) - ^{xV)\{x,t)y. (iA;)^(iA;)^P(x,i)e-*^'^(^'*)x 

X {P^p{k)g^ + P^^ik)gT)\kr''h{x) 

= -R{x,t)^C{x,t)-'^/'-' h{x) (^t^+^-^^L. (3.28) 

For the contribution proportional to , we have retained from the tensor operator 
r{x,t){x'VYr{x,t) only the diagonal contribution ^r{x,t)x'^{A)r{x,t), which is sufficient 
at 1-loop order. For the subtleties associated with the insertion of this operator at the 
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Figure 3.1: The diagrams Di (left) and D2 (right). 
2-loop level cf. [52]. Using the FDT, Eq. (3.12), this can still be simplified to 

|c(.,r"'^(5™+^-^)«L. (3.29) 

Eqs. (3.24) and (3.29) contain all possible divergent terms in the short-distance expansion 
of Eq. (3.20) and all terms which have to be taken into account in 1-loop order. Notably, 
due to causality, no term independent of f appears. 

3.4 The renormalization of the disorder (divergences associated with bilocal 

operators) 

There are also UV-divergent configurations associated to bilocal operators, which renor- 
mahze the disorder, and which are depicted on Fig. 3.1. A dashed fine indicates that the 
end-points of the interaction vertices are approached in space and time. Up to permuta- 
tions of the two interaction vertices, there are two possibilities to order their end-points 
in time, namely Di and D2 on Fig. 3.1. We first calculate Di, starting from 

Uy, t) I {gTPlpik) + gi^p^pik)) \kr'h{x - yy^^<y^')-^^^^')^fp{x, o) x (3.30) 

J k 

X f,{y', t-a)J^ {qtPUp) + ^^^^(p)) W'^Kx' - y')e^^^<y' ^-^^^n{x\ -r) . 

For small x — x', y — y', r and a, with r, a > 0, there is one contribution for the renor- 
malization of the interaction vertex. First, due to causality, r^{y',t — a) and rs{x', — r) 
have to be contracted with a correlator field in order to obtain two response fields at the 
end. Then the short-distance expansion for nearby vertex-operators reads: 

^ikr{y,t)^ipr{y' ,t-a) ^ ^i{k+p)r{y,t) ^kpC{y-y',a) 

This yields up to subleading terms: 

/J^r.(y,^)r>(x,0)e^(^+^)(^(^'*)-'^(-'O)) {gTPlp{k) + g^P^^{k)) {g^Pliv) g^P^^sip)) 

{ik)^{-ik)sR{y - y',a)R{x - x', r)e^f(^(2'-f''^)+^(=^-^''^))|A;|'^-'^|p|"-'^/i(a; - y)h{x' - y') . 

(3.32) 
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In the next step, first k and second p are shifted: 



k-p , 



p 



p + 



(3.33) 



The result is 



/ raiy,t)ff3ix, 0)e' 

Jk 



k{r{y,t)-r{x,0)) 



{q^PUp - k/2) + g^P^pip - k/2)) 



a—d 



(g^Plip + kl2) + g^P};,{p + k/2)) {k/2 - pUk/2 - p)s\p + k/2r'\p - k/2\ 

R{y - y', a)R{x - x', T)e^^V^-p')i(^^y-y' ^'^)+<^i^-^' ^^))h{x - y)h{x' - y')] . (3.34) 

To compute the correction proportional to the disorder vertex, the expression in the 
rectangular brackets is expanded for small k. As the integral has a well-defined limit for 
/c — >■ 0, convergent for d > 2, neither a term of the form k^kf^/k"^, nor a term of the form 
\k\"'~'^ with a < d can be generated. We conclude that in the LR-case only SR-disorder is 
generated. In the remainder of this sections, we therefore focus on short-range correlated 
forces setting a — d. The leading term of the above expansion is then (the p-integral 
being defined in Eq. (2.12)) 



a0 - P "PaPd ) + QhP "^PaPd ) fi'LP'e 



xR{y - y', a)R{x - x' , T)h{x' - y') 



J2^-p'^(C{y-y',a)+C{x-x',T)) 



X 



d 



+ 91^2 ) 9h X 



d 



X ^ {C{y - y', a) + C{x - x' , r))"'^/'"^ C{y - y', a)C{x - x', T)h{x' - y') 



9t 1 



+ 9l-^ ) 9l X 



d 

x^-^£- {C{y - y', a) + C{x - x', r))-'"^' h{{x - y) - {x' - y')) . 



d-2dTda 



(3.35) 



In the last formula we have rearranged the product of /i-functions, using the fact that 
they are always ^-distributions for the cases of interest. Note also that we have used the 
(perturbative) FDT, Eq. (3.12), for the first transformation. 

The second possible way to do the contraction, see D2 on figure 3.1, is performed 
similarly. The leading term is 

~^-9l-A^izi: {C(y - y', a) + C{x - x', r))-'"'-' h((x -y)- {x' - y')) . 



d d — 2dT da 



(3.36) 



Note that there are two other possible contractions, which can be obtained from Di 
and D2 by replacing r and a with —r and —a respectively. Together they add up to 



gL9T 



d-2 



1 - 



)§^§^iC{y-y',a) + C(x-x',r))-'/'^'x 



dj dr da 



xh{{x-y) - {x' -y')) . 



(3.37) 
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This result is remarkable in several respects: First, the contribution to the renormalization 
of disorder from the disorder-disorder contraction is isotropic. We will see that this 
stabilizes the isotropic fixed point. Second, there is no divergent contribution in the 
purely transversal or purely longitudinal case at 1-loop order. (Note however that there 
are finite contributions in the transversal case, see appendix C, such that isotropic disorder 
will be generated at 2-loop order.) 



3.5 Beyond the leading order 

The MOPE-coefficients evaluated in the last two subsections arc sufficient to derive the 
RG-equations at leading order. Subleading orders are analyzed similarly, with the result 
that the only diverging contributions, and thus those which have to be renormalized, are 
proportional to the terms already present in the action. An explicit example is the second 
order term from the contraction of three disorder- vertices towards a single one as 




Another question is, whether the fixed point as given in leading order, is stable towards 
perturbations w.r.t. subleading terms, i.e. those which are irrelevant at the trivial (Gaus- 
sian) fixed point. This is a difficult question, and cannot be studied in its generality here. 
However, we shall show later in section 4.7 that the physically important perturbation by 
self- avoidance is indeed irrelevant at the 1-loop fixed point to be described below. 



3.6 List of the MOPE-coefficients 



We have identified the 1-loop divergences, and shall now compute explicitly the counter- 
terms needed to render the action finite to 1-loop order. To summarize, we first give a list 
of the corrections to each term in the original action in the form of the MOPE-coefficients 
for the diagonal case {gi, — qt). Then the MOPE-coefficients for the transversal and 
longitudinal disorder are expressed with its help. To simplify notations, the fluctuation- 
dissipation theorem (3.12) has been used. Concerning the notation: {A\B) means the 
contribution from the contraction A proportional to B, and is expressed in terms of the 
relative distances. Prom Eq. (3.24) we obtain 





The other local divergences arc determined from Eq. (3.29): 



(3.39) 
(3.40) 
(3.41) 

(3.42) 
(3.43) 
(3.44) 
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and 



2dD dt 



. 



(3.45) 
(3.46) 
(3.47) 



Diagrams proportional to , renormalizing the stiffness, will of course vanish in the 

directed case. 

The bilocal divergences, proportional to the disorder vertex, are 




V\A# •XAA, 



wx* •vw 




dj dr da 





. 



(3.48) 
(3.49) 

(3.50) 

(3.51) 



Remember that in the LR-case, this class of diagrams vanishes. 

3.7 The residues in the directed case 

Let us start by analyzing the directed case, i.e. h{x) — S^{x). Using the MOPE coefficients 
given in the last subsection, we can calculate the diverging part of the diagrams as pole 
in 

£ = 2 - Coa , (3.52) 

where 



2-D 



(3.53) 



and setting a = d for SR. For e > 0, the integrals arc UV- convergent and we put an 
IR-cutoff L in position space and in time. Note that reintroducing A this would yield 
a "physical" time cut-off XL^. The precise procedure is irrelevant at 1-loop order, cf. [53]. 
Note that 

h{x) = S^{x) (3.54) 

always eliminates one space integration. 

Start with an explicit example. Let us define a Feynman-diagram as the integral of the 
MOPE-coefficient over all distances in space and time, bounded by L and respectively 
(for details, see [52]). For the MOPE-coefficient (XJ^I— ) this reads 



L2 



/ diC(0,i)-"/2 

^0 
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Jo 
.2 



P-L' 
e 



where 



a/2 



(3.55) 



(3.56) 



We introduce the foUowing abbreviation for the pole-term, which has to be subtracted in 
the minimal subtraction scheme 



A 



dir ._ 



^AA•vw ) ^P 



(3.57) 



The divergent contribution to the friction coefficient is obtained from 



i/2 



L Jt<L'^ 

t d 



= - di^c(o,i)-«/' + e»(£°) 

Jo a 

= --P r dU-^°"/2 ^ C)(£°) 
a Jo 

^--ph' + 0(e') , 

a e 



(3.58) 



with P introduced in Eq. (3.56). The pole-term, which has to be subtracted in the minimal 
subtraction scheme, is 



Remarkably 



(3.59) 



(3.60) 



This relation can be viewed as a consequence of the (perturbative) FDT, and is responsible 
for the non-renormalization of the temperature in the longitudinal case, where a FDT 
holds in the full theory. 

Since the disorder- vert ex is local in x-space ("statistical tilt invariance" ) , there is no 
contribution proportional to f(— A)r: 



. 



(3.61) 



The renormalization of the disorder interaction is (for a = d only) obtained from 



a/w* •wv 

L2 







X 

dr / da ^ f 1 - 3) l-l- {C{x, t) + C{x, a))-'/'^' 



d-2 
^ dx 



dj dr da 



-] Sd I — (2C(x, 0))-'^/'+' + finite terms , 
dJ Jo X 



(3.62) 
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D 

Figure 3.2: The function I{D). The dots have been obtained numerically, the interpolating line is a 
polynomial fit of fifth degree, later used in the extrapolations. 



resulting in 



where 



^dir 



l-d/2 



Q 



s 



D 



(3.63) 



(3.64) 



^{2-D)Sd) 

with Sd = 27r^/^/r[D/2]. This terminates the evaluation of the Feynman-diagrams in 
the directed case. 



3.8 The residues in the isotropic (non-directed) case 

In the isotropic, i.e. non-directed case with 

h{x) = 1 , 

the dimensional regularization parameter is 

£ = 2 + D - Coo . 



(3.65) 



(3.66) 



We now parallel the calculation of diagrams from the last subsection. The term propor- 
tional to (see Eq. (3.39)) is 



L JQ JO 



d^x I dt 
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dx 




r 






Jo 


X 






roo 










£ . 


10 



X 



D 



Jo 



(3.67) 



where in addition a has be set to ac{D), defined by e{D, ac) = 0. Here Cc = 2(2 + Z))/(2 — 
D). The residue is 



= Sd dtC(l,t)-'^=(^)/2 
Jo 



(3.68) 



where the expression for the free correlator (obtained e.g. via time integration of the free 
response function) reads 



C(l,i) 



r 



+ 



2(2 - D)t:2 An 



r -1 + 



D 1 
2"' 4t 



1 



D 



r(^)+(4i)'^e- « 



2(2-D)7r2 

where we have abbreviated the incomplete F-function as 

r(a,6) 



^ 2 ' 4t 



(3.69) 



(3.70) 



Wc could not find a closed analytical form for the integral I{D). Using the approxi- 
mation for the correlator 



C{l,t) 



, 2 



(2-L>)(47r)^/2 

which works astonishingly well numerically, we obtain for the integral I{D) 
/(.) = ..(^^^«)"(|)^""%~ 



(3.71) 



^SD{SD{2-D)f''' X (2^-'r 



D) 

DJ 



(3.72) 



where a must be set to a = ac{D). If the approximation in Eq. (3.71) was exact, I{D) 
would equal 1. This is not the case, and the value for I{D) obtained from numerical 
integration is plotted on figure 3.2. It is a slowly varying function though, well behaved 
at D = and D = 2. Note also that in Eq. (3.72) a factor of Sd {Sd{2 - D)f^'^ is 
factorized which is common to all other diagrams and will therefore cancel at the end. 
Proceeding equivalently, the diagram correcting the friction coefficient is calculated as 



a Jo 



(3.73) 
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Like in the directed case, we have 

XL 



-d 



(3.74) 



This relation can be used to simphfy the RG-calculations. As in the directed case, it is 
a reflection of the FDT, vahd in perturbation theory, and in the longitudinal case in the 
full theory. 

In contrast to the directed case, the stiffness of the membrane is renormalized through 



dt 



2dDJo Jo dt ^ ' ' 

= -^^^ t C{x, 0)-'^/' + finite 

2dD Jo X 



2dD 



ljzSn{{2 - D)SoT'^— + finite 



(3.75) 



The residue is therefore given by 



■■={£L-~^)=-^^So((2-D)Snf' . 



i - 'u UD - (^-^^ 

Finally, the diagram correcting the disorder vertex is (for a = d only) evaluated as follows: 



L 

L2 



d- 


2 


4 




d- 


2 


4 




d- 


2 


4 




d- 


2 



^--)sl — {C{x, 0) + C{y, 0))-'^/'+' + finite terms 
dj Jo X Jo y 

1 \ poo Aq* I T ^ 

1-^)SI — (C(x, 0) + C(l, 0))-"^/'+' — + finite terms 

1-^)SI ((2 - D)Snf'-' r -^"^ f^'"'' + l)"'^'^' - + finite terms 
dJ Jo X ^ ' e 

(3.77) 



2-D 



+ finite terms 



The residue is given by the expression 

4 



e d- 



1 

1 - - 



2V d){2-DyT{^) 



X Sd ((2 - D)Sd) 



d/2 



(3.78) 

The factor involving the F-functions is familiar from the renormalization of the interaction- 
vertex for self-avoiding membranes. 
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3.9 The particle limit 



In this section, we demonstrate that our model reproduces, both in the directed and 
in the isotropic case, the known results for the particle upon taking the limit D ^ 0. 
Some subtleties have to be taken into account. Study first the directed case. Then 
unambiguously 

= 2 . (3.79) 



I e 



lim ( (J 

The other independent and non- vanishing digram is in the limit of small D and e 

2D 



\ dir 

/ £ 



d-2 



(3.80) 



Note that it was calculated under the tacit assumption that e ^ is taken first. In order 
to perform the limit -D ^ 0, we insert d — 2 — 2^^ into Eq. (3.80) and then perform the 
limit £ — > 0. We obtain 



D{2-D) ^ 
D-e ' 



2-D ^2 



(3.81) 



One easily verifies that these results coincide with a direct calculation for the particle, 
obtained by using the single particle propagator and response functions and by omitting 
the space integration from the very beginning. The isotropic case reduces to the particle 
case along the same lines. First 



lim 



n 



V 9 

I £ 

\ iso 

= 



(3.82) 



The last diagram is 



iso AD 
d-2 



(3.83) 



Note the apparent difference in factor of 2 to Eq. (3.80). Inserting now d — 2 — 2 
yields 



2D- 



2-D 



-VXA* •WV 



2D(2-D) ^ 



2D-e 

and thus the same final result, cmp. Eq. (3.81). 

Our calculations, at least at this order, thus always include the particle case. 



(3.84) 



3.10 The flow-equations 

We now introduce renormalized quantities. Going back to the original model, Eq. (3.2), 
we take care to reinsert the factors of A, which were dropped in intermediate expressions. 
Define 

ro = ^^Z r 

fo = /i-f (3.85) 
^0 = Zx \ . 
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Explicitly, the counter-terms are as follows. First, introducing a counter-term for 
yields 



9t 



'WXM- ) — |— 



\ 



9t 



9l 



'ZZ^ 1 + 
= 1 + 

The counter-term for A is given by 
In the isotropic (non-directed) case, there is also a counter-term for A namely 

\fz^z^ = 1 + / Qt 

= 1+ ' 



£ e 



e 



These equations are solved for Z, Z and Z^: 



Z 



1 + 
1 + 

1 + 

+ 



) — 

/s \d s 



) + 



1 - 



e 

£ 

1\ 9t 



dj e 



l9L 

d e 



+ 1- 



/e/ e 



(3.86) 



(3.87) 



(3.88) 



(3.89) 



(3.90) 



(3.91) 



It is more complicated to introduce renormalized couplings gi, and qt- Set 

(pL^^L + PI09t) AV -a ^ (P>° + P>°) A^ . 

This equation is to be understood such that quantities on the l.h.s. are renormalized, and 
those on the r.h.s. are bare. Solving the system of equations yields 



(3.92) 



= i.^Z-^Z'^I^Z-' [g^ - \ 
g% = ifZl^Z^I^Z-' [g^ - 




'W^ •vw 



9l9t\ 



£ e 



j + 0{g') . (3.94) 



We are now in a position, to define the so-called /3-functions, quantifying the fiow of the 
renormalized theory upon a variation of the renormalization scale, through 



9l 



9t 



(3.95) 
(3.96) 
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The /9-function for the longitudinal disorder coupling is 

d 



/5l(c/l,^t) = 1^ 



djjL 



9h 



(3.97) 



W\* •VIA/ 



£ e 



-SQL + 



a + 2 



a + 2 
2d 



+ 



a - 2 



1 



WV* •wv 



Using the relation 



(3.98) 



(3.99) 



valid both in the directed, see Eq. (3.60), and isotropic, see Eq. (3.74), case, Eq. (3.97), 
results in the simplification 



PhigugT) = -SQL + 



a-2 



£1 



/. 2 



-WV* 0V\Aj 



+ 0{g') . 

By an analogous calculation, we obtain for the /3-function of the transversal disorder 
d 



(3.100) 



dfjL 



9t 



(3.101) 



'W^ •vw 



9l9'T..-2s 



£ e 



= -£9t + 



a + 2 



a + 2 
2d 



+ 1- 



a-2 



1\ a + 2 
d 



/ e 

■w^*w,. ) ] + 0{g 



e 2 

' ^ 2 I rn/„3\ 



I £ 



(3.102) 



With the simplification given by Eq. (3.99) the final result is 



/5t(5'l,5't) = -egT: 



a-2 



— ^ ghgT 



^,^_na + 2 



9't + 0{g') . 



(3.103) 
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In view of a possible stable fixed point for gi^ = g^^, it is interesting also to study the 
flow-equation for the difference of the both couplings gh — Qt- 



a-2 



H-)5'L(fi'L - Qt) 



— )9t{9l - Qt) + 0{g^) . (3.104) 



3.11 Other RG-functions and critical exponents 

Several critical exponents can be identified. Taking into account scaling relations among 
them, there are three independent exponents in the directed case, and four in the isotropic 
case. The first is the full scaling-dimension of the membrane, given by the fixed point 
value C* of 



C(^l,c/t) 



2-D Id 
Li—- InZ 

2 2^a//„ 



2-^ Ql 
2 2 



+ 0{g') . 



(3.105) 



We again used Eq. (3.99) to simplify the result. The dynamical exponent z is 

d 



InZ, 



( d d ^ 

2 - ( /3L(^L,^T)^lnZA + /3T(^L,^T)^lnZA 



2 + 



The (full) scaling dimension of the response field is 



+ 0{g') 



(3.106) 



C{9u 9t) 



2+D I d ~ 

11-— InZ 

2 2^d^ „ 

\ (/9L(5'L,5'T)^ln^ + /3T(5'L,fl'T)Tr-ln^ 



2 2 
2 + £> 5l 



/ e 



We used Eq. (3.99) to simphfy the result. 



w*w>,y -)- Oi^g ) . 



(3.107) 
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In the isotropic (non-directed) case, there is still another independent exponent, which 
is called j3 (not to be confounded with the /3-functions describing the flow of the coupling- 
constants) and which measures the anomalous dimension of the stiffness of the membrane. 

= ^l(X1|-^)^ + 0(^^). (3.108) 

With the help of the above introduced exponents, we can quantify the non-linearity 
in the response of the drift velocity to a uniform applied force. Adding a small constant 
term / to the random force in Eq. (2.1) results into a change in the dynamical functional 
of the form 

A /' f(x,t)/. (3.109) 

Jx,t 

Since this term does not require a renormalization, we deduce the exact exponent identity 
for the full scaling dimension [/] of / 

[/] + C + ^ + i?-0. (3.110) 

From the definition of /3 in Eq. (3.108), we have a second identity, namely 

/3 = C + C + ^ + L»-2. (3.111) 

Combining these two relations, we obtain 

[/] = C-/3-2. (3.112) 

Since the drift velocity has dimension 

H = 



(3.113) 



we deduce the scaling relation 



z — C 

f with = . (3.114) 

Since, as we will see, always z > 2 and /3 < 0, is larger than 1 indicating trapping of 
the membrane by the random flow. 

Let us also mention another possible parameterization of our model, where one does 
not introduce A and renormalization factors for A and the field r, but associates a friction 
coefficient 77 with r, entering into the dynamical functional as 

T] / r{x,t)r{x,t) (3.115) 

Jx,t 

and a temperature T associated to disorder, entering into the dynamical functional as 

riT [ f-^{x,t) . (3.116) 
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The corresponding renormalization factors then read 




(3.117) 
(3.118) 



It is especially enlightening to consider the renormalization of temperature, which in 
conserved systems (i.e. those with longitudinal disorder) is preserved by the fluctuation- 
dissipation theorem. Stated differently, temperature appears as a parameter in its static 
partition function, and dynamics is constructed such that via detailed balance the static 
partition function is reproduced. Only non-conserved forces can give rise to a change in 
temperature. This is read off from the explicit form of at 1-loop order: 



Oh 

e 




(3.119) 



Since the term proportional to gi^ vanishes as found in Eqs. (3.60) and (3.74), anomalous 
corrections to the temperature are proportional to and present whenever g-^ 7^ 0. 
Stated differently: In the presence of non-potential forces, the fluctuation-dissipation 
theorem (FDT) is violated. 



4 Results and discussion 

In this section we analyze the general renormalization group flow given by the two f3- 
functions in Eqs. (3.100) and (3.103) where the coefficients have been computed in section 
3.7. We identify the fixed points and compute the critical exponents at these fixed points. 
Before doing so it is useful to recapitulate the results obtained so far in the previous 
sections. 



4.1 Preliminaries 

To present conveniently the analysis below we introduce the following notation for the 
three independent coefficients (two in the directed case) as 




The explicit expressions for these coefficients are given above: for A see Eqs. (3.57) and 
(3.68), for B see Eqs. (3.63) and (3.78) and for C (only non-zero for isotropic manifolds). 
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see Eq. (3.76). They are non negative. In terms of these three coefficients, the /9-functions 
of the two coupling constants read: 



(^l{9u 9t) = -SQL + ^ ((l - (a + 2)^ - qi^qt - 
/3|(^L, ^t) = -SQT 9T9L + (l - ^ 



a + 2 



(4.4) 
(4.5) 



where a must be replaced by d for short-range disorder (the crossover between SR and LR 
disorder is also studied below). Similarly, the expressions of the exponents are to lowest 
order (from section 3.11 above): 



z{gL,gT) ^2 + -{A-C)gi^ 



(4.6) 
(4.7) 

(4.8) 



where we have used (3.74). Note that in the end the exponents will depend on the 
amplitudes only through their ratios BjA and CjA which are universal, in contrast to 
the amplitudes themselves, which depend on the specific normalizations of the model. 

The coefficients of the RG equations (4.4)-(4.8), including A^ B and C, depend a priori 
on a, d and D, and their full dependent expressions are given in the previous sections 
(where SR and LR disorder are treated in a single calculation). However, since we are 
expanding to first order around the line of critical dimension dc{D) (see Eq. (2.22)), we 
can simplify all these coefficients as follows: For short-range disorder one can substitute 
a — d — dc{D) in all coefficients. For long-range disorder we can substitute a — ac{D) 
everywhere, but one must keep d as an independent variable. This is the reason, why we 
defined A and C as d independent, while i3 = for LR disorder. The final result for the 
coefficients used below, expressed only as a function of D is, for directed manifolds: 



A'^" = 2 



(2-L>)(47r)^/2 



2 

2-D 



2 + D 



D 



r(f) 



2 

2-D 



(4.9) 



7 dir 



. 



For polymers {D = 1) one finds A — 27r, B — 3 and the universal ratio A/B — 27r/3. 
For isotropic manifolds one has: 



^iso 

^iso 



2-D 



(2 - D)22+-D \Dj 

2D{2-DY{2 + D) r(^) 
1 



G 



D 



(4.10) 



2D{2-D) 



G 
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where we denote by 

the factor common to the three amphtudcs (which therefore drops in the final results). 
I{D) is plotted in Fig. (3.2) and defined as 

(/") \ — ~ — /'+00 2+D 

" 2^+^/^ dt[r(f) + (4t)^e-i.-r(f,i,)]--- . (4.12) 

It takes the values /(D = 0) = 1, i{D = 1) = 1.82202 and i{D = 2) = 3. For the 
coefficients ior D ^ 1 we find A = 1.79351Gi, B = |Gi and C = ^Gi. 

We now specify to each of the four cases of interest and compute the critical exponents. 

4.2 Directed manifolds with short-range disorder 

The ffow for this case is represented on Fig. 4.1. The expansion parameter is £ = 2 — (2 — 
D)d/2. The fixed points arc found to be: 

(1) Gaussian fixed point 

The Gaussian fixed point at (^l = 5't = is completely unstable for d < (i^"'(i5) 



2-D- 



(2) Potential disorder 

The fine = is found to be preserved by RG and we find a flow towards strong 
coupling. This case corresponds to the dynamics of a directed manifold in a long-range 

correlated random potential (short range correlated force). It has been well studied and 
it is indeed known that for d < dc{D) (defined above) and D < 2, or in any d for 
2 < D < 4, the physics is controlled by a strong disorder (zero temperature) fixed point 
[1,6], inaccessible by the present method (it is accessible only via a, D — A — e calculation). 

Further results are available for the particle {D — 0, dc — 2). There it is known 
rigorously that for = 1 (Sinai model) ultra-slow logarithmic diffusion r ~ (Int)^ occurs 
(z/ = 0) and it is believed (but not proven) that similar ultra-slow diffusion also occurs 
for manifolds D > (see e.g. [12]). At the critical dimension, for the particle dc{D — 
0) = 2, 1- and 2-loop calculations [39,40,42] have shown that diffusion is anomalous with 
a continuously varying exponent iy{gi,)- It was further argued that the /3 function vanishes 
to all orders, and that the one loop result for i^igh) holds exactly, for SR disorder as well as 
for LR disorder [40]. While unchallenged for SR disorder (i.e. at d = 2), the exactness of 
the one loop exponent was recently questioned for LR disorder in [76]. To our knowledge 
such issues have not been addressed for manifolds. 

(3) Isotropic disorder fixed point 

Noting from Eq. (3.104) that the line qt — Ql is preserved by the ffow at 1-loop order, 
we find an isotropic fixed point at 

+ (4.13) 

where the coefficients are given in Eq. (4.9). Since at the critical dimension the denom- 
inator is positive for all values of < < 2, ^f* is positive and thus physical. This 
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Figure 4.1: RG flow diagram for SR disorder. The physics is controlled by the fixed point /* at (?t = dh- 



fixed point is completely attractive, and its domain of attraction covers all perturbative 
situations except the potential case with — 0. As we show below, it also controls the 
line (^L = (except for D = 0). 

Let us give the values of the roughness exponent (^* and of the dynamical exponent 
z*, defined in Eqs. (2.13) and (2.15), at the fixed point: 



B 



2-D 



e 



R 



2 d+2-R 
d 

(1-1)^-^-2 



4- (2^-r(i;/2)) 



d/2 



(4.14) 
(4.15) 



where e — 2 — ^-^d. Again one can substitute d — dc{D) in the coefficient of e. Expressed 
as a function of D one has R — ^{2^~^V{D /2))'''^/^'^''^\ The dynamical exponent is: 



2 + 



2£ 



{d-l){d + 2-R) 



(4.16) 



The exponent (3 equals as a consequence of the statistical tilt symmetry (see Eq. (3.61)). 
The exponents v and can then be found using 







2-C 



(4.17) 



In particular, for the polymer D = 1, we find that below rfc = 4 the random fiow is 
relevant and roughness and dynamical exponents take the anomalous dimension 



^ 2\ 3-7r-i 



= 2 1 



6(3-7r-i)^ 
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(4.18) 
(4.19) 
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Figure 4.2: Results for directed polymers and membranes, isotropic fixed point, SR-disorder. Values given 
for the exponent C (obtained by extrapolating Cd) arc rather faithful. Corrections to z arc small, but no 
plateau could be found to optimize the extrapolations. The given values are therefore only intended to 
give a rough idea. The same is true for extrapolations for v. We have reported the values of v obtained 
by setting Dq = D. 



where e = |(4 — rf). Interestingly we find that to first order in e, ( > (p. 

The most naive fixed D expansion gives in d = 3 (setting £ = 1 in the above) that 
C = 0.686, z = 2.12, 1/ = 0.323, = 1.094. In d = 2 it gives (setting £ = 2) that C = 0.873, 
z — 2.249, u — 0.388, cf) — 1.22. These estimates can be improved since the expansion can 
be carried from any point of the line d = dc{D). It is possible to optimize the numerical 
estimation of the above critical exponents over the expansion point Dq. Details of the 
general optimization method are described in [52] and applied to the present case in 
Appendix B. The results are listed in table 4.2 and are purely indicative: we have even 
refrained from giving an error bar. 

(4) Transversal disorder 

There is an apparent fixed point at 

where the coefficients are given in Eq. (4.9). Interestingly, at this fixed point we recover 
exactly the Flory-result at 1-loop order: 

z^2 + 0{e^) . (4.22) 

In the case of the particle [D = 0) the result v — v-p was found to be true to all orders 
(see e.g. [42,40] and discussion below). 

However, for D > we show (see Appendix C) that finite non-transversal terms are 
generated in perturbation theory. Since the above fixed point is unstable to adding a small 
longitudinal disorder, > 0, these terms will drive the system towards the isotropic fixed 
point discussed above. 
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Let us close this section by noting that we recover the known results for the particle 
D = 0, obtained via different RG techniques in [36,39,40,42]. That this should be so was 
discussed in section 3.9 but here we discuss it explicitly for the exponents. First we note 
that the flow diagram of Fig. 4.1 has the same structure for all D as the one found in 
[39] for the particle with an isotropic disorder fixed point and the transversal fixed point. 
Second, expressing in all our results the coefficient of e as a function of D alone (using 
d — dc{D)) yields in the limit D ^ with e = 2 — d, that at the isotropic fixed point 

V — C,/e — l/2 + 0{e^) and (p = 1 + e (note that one also finds C = 1 + £/2, z = 2 + e, but 
only V and have a physical meaning for D — Q). Similarly at the transversal disorder 
fixed point (which has a meaning only for the particle), one recovers the known result 

V — and 0=1. 

4.3 Directed manifold with long-range disorder 

Long range correlated disorder (a < d) can be studied in an expansion near = 4/ (2— Z)), 
in the small parameter 5 = 2 — (2 — D)| (with d and D fixed). We have shown in section 
3.4 that in this case the vertex renormalization vanishes 

^ = . (4.23) 

Thus the renormalization is driven by the corrections to temperature and friction. The 
/3-functions in Eqs. (3.100) and (3.103) lead to the flow diagram presented on Fig. 4.3. 
Since there is no vertex renormalization, both couplings have the same renormalization 
factor, the ratio k — Ql/ Qt is preserved by RG, and the flow is along straight lines. One 




Figure 4.3: RG flow for directed manifolds in LR disorder. The flow is along straight lines and preserves 
the line L of fixed points. 

flnds a line L of flxed points, which is vertical to flrst order in e. This feature, already 
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noted in the particle case, D — [42,40], thus extends to manifolds with D > 0. Here the 
fixed line is given by: 

2r/ 

g* ^ S- -j^ —r (4.24) 

{d-l){a + 2)A ^ ' 

For LR correlations the result for C* at any point on the fixed line can be derived to 

all orders in perturbation theory. First, there is no counter-term to fAr to any order 
in perturbation theory due to the fact that the interaction vertex is a difi^erencc of the 
r-fields at the same point ("statistical tilt symmetry''), resulting in the scaling relation 

C* + C + -2* = 2 - D , (4.25) 

where C, is the full scahng dimension of the response field. Second, as long as there is no 
renormalization of the vertex, 

L> + 2C + 2z* - aC* = (4.26) 

is valid at the fixed point. Combining these two relations gives^ 

A- D 

C = Cf^ ^ ■ (4.27) 

^ ^ 2+a ^ ^ 

Thus the Flory estimate for ( is exact in the LR case. Note also that ( in that case does 
not depend on the type of disorder (transversal or potential). For a polymer D = 1 we 
obtain that ( = 3/{2 + a) for a < = 4. This leads to the interesting prediction that 
the polymer will become overstretched (^ > 1) when the range of the force correlation 
increases beyond a < 1. 

The dynamical exponent z, however, and consequently and are non-trivial, and 
depend continuously on the ratio K — g^/gi,. They are given by: 

z*^2(l + K- ^- — -) . (4.28) 

V i2 + a)id-l)J ^ ' 

where 5 = 2 — (2 — Z})|. This result is valid to first order in 5 and one may thus replace a 
by Oc = 4/(2 — D) in the above formula. It shows that the more potential the dynamics 
the more glassy it is (the higher z). 

We can now compare these results with the results obtained in [44]. There the ex- 
ponents were obtained exactly for d ^ oo and arbitrary a. A Hartree approximation 
was obtained for finite d. The notations of [44] were different. To translate to the 
present notations one must set in [44] 7 —>■ a/2, g2 —>■ Qi — §2 (l~^)fl'T and 
£ = 2/(2 — D)— 7— > 23d- With this correspondence the results of [44] are as follows. 
First C = Cf, both for d = 00 and for finite d within the Hartree approximation. This is 
in agreement with the present result, z in [44] however is a highly non trivial function 
of a, D (and of d in the Hartree approximation) implicitly given by equation (9) in [44]. 
Remarkably, for D < 2, expanding z near a = ac, which yields formula (12) of [44], one 

^This identity also holds at 1-loop order at the apparent SR transversal fixed point discussed in Eq. 
(4.21), since vertex-corrections are proportional to ^lS'Ti see Eq. (3.37). 
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finds to first order in e exactly the same result as here in Eq. (4.28). Of course tlie two 
also coincide, as tlicy sliould, for d = oo. Tlius we find tliat in tlie long-range directed 
case, the Hartree approximation in [44] gives the correct result for z to first order, and is 
exact for (. 

Let us also give the results for the exponents u and (p. One finds: 



2-D 



• (k*-k){2-D) 
^ 2(4-L>)(d-l) 



(4.29) 



and thus we find that the single monomer dynamics is either subdiffusive if k, = gi^/ > 
K* — 2{d — 1) / {2 — D) or superdiffusive if k < k*, while is always larger than 1 (except 
when the flow is purely transversal). One can also check that our formulae agree, in the 
limit D ^ with the known formulae [40-42] for the particle (with = 2 and 5 = 2 — a). 

To conclude note that in the case k = 0, i.e. in presence of LR transversal disorder 
alone, the mechanism of generation of potential disorder and barriers is different from the 
SR-case. LR non transversal disorder is not generated, only SR non transversal disorder 
is. (See also the section below, where the crossover from LR to SR is discussed). The final 
result is that in that case the exponents are given correctly too by the above formulae 
setting K = 0. 



4.4 Isotropic (non-directed) manifold with short-range disorder 

For an isotropic manifold, we find from the RG equations (3.100) and (3.103) and the 
expressions given in section 3.8 for the coefficients that the RG fiow is qualitatively similar 
to the one in the directed case depicted on Fig. 4.1, with the following fixed points: 

(1) Gaussian fixed point 

The Gaussian fixed point at g'L = 5't = is completely unstable for £ = 2 + D — |(2 — 
D) > 0. 

(2) Potential disorder 

The line g^ = is preserved by RG and again, we find a fiow towards strong coupling, 
which is faster than in the directed there is an additional term. This problem de- 

scribes the dynamics of an isotropic manifold in a long-range correlated random potential 
(short range correlated force) . The statics of this problem has been much studied and it 
is indeed expected that this problem is described by strong disorder. 

In the absence of self-avoidance some features of the problem are relatively simple to 
analyze. For instance can a polymer (D = 1) in a rf-dimensional random potential also 
be seen as a directed polymer in a d -|- 1 dimensional random potential (the additional 
coordinate being the internal dimension) which is constant (i.e. correlated) in the internal 
dimension, a problem extensively reinvestigated recently in relation with flux lines in 
superconductors in presence of columnar disorder [81,77,78,12]. It is also related to the 
problem of quantum localization of a particle in a dimensional random potential (the 
internal dimension being the quantum imaginary time). It is known, and rather obvious 
in relation to these problems, that if the polymer is free to move, it is localized in the 



36 



deepest localized state and has a finite extent. Properties then strongly depend on details, 
such as the system size, the tail of the disorder distribution (Lifschitz tails), and the short 
scale cutoff [79,80]. 

On the other hand, the problem of many interacting directed polymers is well defined 
and identical to the problem of localization of bosons in a random potential. These 
properties should remain true for membranes as well and in general we expect that in the 
absence of self-avoidance, the membrane collapses = 0). 

For a self-avoiding polymer the problem is better defined. Self- avoidance is then 
relevant (see section 4.7) and properties will depend on whether one end of the polymer 
is held fixed or is allowed to move freely. It was studied previously with short-range 
potential disorder [55]. There, if the polymer is free to move one expects again [60] a 
collapsed state, whereas if one end is fixed there may be another strong disorder fixed 
point [57] (see section 4.7). In all cases it is natural to expect non-trivial glassy dynamics. 

We will not say more here about the potential case (qt — 0) which is beyond the 
reach of our perturbative method, but let us stress that it is all the more interesting to 
investigate what happens to localization when transverse disorder is added. We therefore 
now turn to the case g-r > 0. 

(3) Isotropic disorder fixed point 

Noting from Eq. (3.104) that the line qt = gj, is preserved by the fiow (to this order), 
we find an isotropic fixed point at 

2ed 

where the coefficients are given in Eq. (4.10). We have checked numerically, that the 
denominator of Eq. (4.31) is always positive, which is necessary for this fixed point to 
be stable and to be in the physical domain. We have also checked numerically that this 
fixed point is completely attractive, and its domain of attraction covers all perturbative 
situations except the potential case g^ = 0. As in the directed case, it also controls the 
fine gL — (except for D — 0) and describes the large scale behavior of an isotropic 
manifold in a random short-range force fiow, see appendix C. 

The critical exponents at this fixed point, (* and z*, defined in Eqs. (2.13) and (2.15), 
are with the same diagrams, given in Eq. (4.10) 

Note the difference with the directed case. The first coefficient in the numerator. A, 
is positive as before, since it arises from the upward corrections to the temperature. 
However, since there is no tilt symmetry any more, the elasticity is also renormalized 
upwards (the polymer tends to shrink to take advantage of favorable regions). This 
produces the second coefficient — C, which is negative. The competition between the two 
opposite effects however gives a positive sum and the membrane is stretched. Also note 
that the e-correction vanishes like 1/d for D ^ 2. 

A similar formula is valid for the dynamic exponent z: 
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Figure 4.4: Results for isotropic polymers and membranes at the isotropic fixed point, SR-disorder. 
Results for C arc obtained from the extrapolation of (d, and are the most faithful. One observes a nice 
plateau in extrapolations for v, but no corrections to the Flory approximation can be deduced. The 
exponents z and ^ do not allow for direct extrapolations, but are always corrected upwards. Results for 
(3 are significant for polymers. For membranes only a bound seems to emerge from extrapolations. 



with the same coefficients. The ^-correction is always positive, but vanishes like for 

Since there is no tilt symmetry any more, the elasticity is also renormalized upwards 
and gives rise to a non-trivial exponent P* (with the coefficients again given in Eq. (4.10)): 

The other exponents can be obtained as u = (/z and = {z — C)/{2 — ( + j3). One 
notes that in the limit D — > one recovers again correctly the results for the particle for 
V and 0. 

For a polymer, D — 1, the disorder becomes relevant below d — & and setting D — 1, 
we find that the above results yield: 

= 0.5 + 0.130792£ , z* = 2 + 0.03996£ , /?* = -0.015446£ (4.35) 

with e = |(6 — (i), as well as z/ = 1/4 + 0. 060401^ and 0= 1 + 0.0369373^. The most naive 
extrapolation to o? = 3 would be (setting e = 1.5) that ( = 0.70, z = 2.06 and (3 = —0.023 
and in = 2 that C = 0.76, z = 2.08 and P = -0.031. 

One can try to obtain more reliable estimates for these critical exponents from expres- 
sions (4.32), (4.33) and (4.34) for polymers {D = 1) and membranes {D = 2) in three or 
two dimensions by optimizing on the expansion point. This is a tedious task since e is 
rather big. The numerical values obtained by the methods of [52] are not very precise, as 
we could not find a combination of the exponents and D or d, which in suitable extrapo- 
lation variables build up a nice plateau. Different extrapolation schemes yielded strongly 
varying results. Some indicative values obtained by this methods are summarized in Fig. 
4.4. 

Since ( seems to increase rapidly as d decreases, an interesting question is whether 
there is a dimension di below which the polymer will be fully stretched {( = 1). The 
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result of Fig. 4.4 seems to indicate that di could be around di = 2. Our calculations are 
not precise enough to decide on whether the polymer is already over-stretched or not in 
d — 2 but that would be an interesting point to try to answer by other methods or by 
numerical simulations (in d — 2). 

(4) Transversal disorder fixed point 

The transversal fixed point (g'L = 0) is at 



where the diagram is given in Eq. (4.10). It is unstable towards perturbations of gi^ see 
appendix C. For the critical exponents (* and z*, we recover exactly the Flory- result at 
1-loop order: 

C* = ^ + dT2 + ^^"'^ ^^-^^^ 

z* ^2 + 0(£^) . (4.38) 

As was discussed for the directed case, this result was also found to be true for the 
particle {D — 0) to all orders, but for D > as we show (see Appendix C) finite non- 
transversal terms are generated in perturbation theory which will drive the system towards 
the isotropic fixed point discussed above. 

4.5 Isotropic (non-directed) manifold with long-range disorder 

For isotropic manifolds, LR correlated disorder becomes relevant for a < ac{D) = 2(2 + 
D)/{2 — D). For fixed dimension d > ac{D), it can be studied in an expansion in S = 
2 + D — ^^~^-* a. As for directed manifolds, the vertex renormalization vanishes. 



B 



— ©wv 



. (4.39) 



Thus the renormalization is driven by the corrections to temperature, friction and elas- 
ticity. The /9-functions (4.4) and (4.5) with the coefficients given in Eq. (4.10), lead to 
the fiow diagram represented on Fig. 4.5. Since there is no vertex renormalization both 
couplings have the same renormalization factor and thus the ratio k = Qh/gT is again pre- 
served by RG. As in the case of directed manifolds, one finds that the flow is along straight 
lines and that there is a line L of fixed points (see Fig. 4.5). Its equation parameterized 
by the bare ratio of coupling constants k — 5'l(0)/5't(0) is: 

^ _ 2dS 
" {d-l){a + 2)A- K{a-2)C ^^'^^^ 

and = '^S'tj where a should be replaced by a = ac{D) = 2(2 + D)/{2 — D). 

A novel and interesting feature is the existence of a separatrix, as depicted on Fig. 4.5: 
When starting with mostly potential disorder, i.e. a large enough ratio: 

K=f>K.:^{d-l)^^, (4.41) 
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Figure 4.5: RG flow diagram for isotropic manifolds in LR disorder. The flow is along straight lines. 
There is a line of fixed points L and an apparent separatrix. C ~ Co changes sign upon increasing sl/sT) 
suggesting progressive localization. 



there is no perturbative fixed point and the system, within the present expansion, always 
flows to strong couphng (for the coefficients see Eq. (4.10)). For polymers D = 1, since 
we expand around a = = 6 this value is rather large Kc{d) = 7.17404((i — 1). This 
property suggests the existence of two phases, or two regimes, one described by the above 
fixed line, the other by a fixed point, which is non-perturbative in 5 (at fixed d). It 
could also be that while the perturbative one has a continuously varying exponent ( the 
non-pcrturbative corresponds to complete localization C = 0. It would be interesting to 
check this scenario by non-perturbative methods as e.g. in [44]. A study of the crossover 
between SR and LR in the next section will shed some more light on the nature of the 
change of regimes observed here. 

In the perturbative (non-localized) phase k < Kc the exponents as a function oi k — ^ 
read: 

^ 2 {d-l){a + 2)A-{a-2)KC ^' ' 

z* ^2 + 26k— "^'f , — (4.43) 

{d-\){a^2)A-{a-2)KC ^ ' 

where a should be replaced by a = a^iD) = 2(2 -D)/(2 - D) and 6 ^2^D - ^a. 
Again the diagrams are given in Eq. (4.10). 

These results contain interesting novel features compared to the case of directed man- 
ifolds and can be studied as follows. Let us denote the universal ratio R = A/C. We can 
express all exponents as follows to first order in 5 — 2 + D — {2 ~ D)^: 

2 — D 6 Kn — n 

c = \ — 

2 Oj 2 tvp t\j 
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25 niR-l) 

Qj 2 K/(i^ AC 

,5=-^^ (4.44) 
a — 2 Kc — K 

^5 Rk 
^='+ („_2)(2 + 0)K-«) 

with = (rf - and = 2{d - l)Rl{D + (2 - and one should set a = adD) = 

2(2 + D)/(2 — D), i.e. a = Oc = 6 for the polymer. 

In the case of the polymer {D = 1) we find that R = 3.58702, and k,c = 2R{d — 1) = 
7.17404((i — 1) varies in the range 35.8702 < Kc < +oo for the range of dimensions 
ttc — 6 < d < +00 where the result holds. Kq — R{d — 1) = Kc/2 = 3.58702(d — 1) and 
K* = 1.56399(ci — 1) which varies between 7.81993 < k* < +oo as d increases from d — 6 
to d = +00. Note that one has k,* < kq < Kc- 

The above given lowest order results for the exponents in the expansion in d = 2 + D — 
^-^a show some new features as compared to directed manifolds which result from the 
competition between the effect of temperature increase (which tends to stretch) and the 
effect of increase in elastic stiffness (which tends to localize). First we find that for k < Kq 
(with Kq = Kc/2 for polymers) the manifold is stretched ( > compared to manifolds 
in the absence of disorder (as was always the case for directed manifolds) while for large 
potential disorder k. > kq the manifold is more localized C < Co- The effect of potential 
disorder, which acts to localize isotropic polymers as was discussed above thus becomes 
predominant for k > kq (to lowest order in 5). The single monomer diffusion exponent 
r t" also exibits a change from superdiffusive (z/ > 1/4 for polymers) to subdiffusive 
(z/ < 1/4) behavior when k becomes larger than k*. Finally note that in all cases k > 
one has z > 2 and (p > 1 which indicates glassy dynamics and trapping by the flow. 



4.6 Crossover from short-range to long-range correlated disorder 

It is interesting to discuss the crossover from short to long-range disorder. It can be studied 
analytically around the point a — d — dc{D) in a double expansion. In the case of the 
particle, this has been done by several authors (see [42] and references therein) in a double 
expansion near a = d = 2. A complicated scenario was found with several sectors. We will 
recover their result and flnd novel ones for manifolds. For polymers we will expand near 
d = a = 4 for directed polymers and d = a = 6 foi isotropic ones. Typically one expects 
a line ac{d) where the SR exponents cease to hold, as represented on figure 4.6. It is thus 
important to study the crossover in detail in view of numerical simulations since in the 
physical dimensions d = 2,3 there will be a wide range of correlation decay exponents a 
accessible in practice, where the competition between long-range and short-range disorder 
obtained here can be tested. 
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Figure 4.6: Schematic sketch of the crossover in the d, a plane between LR and SR regimes for a directed 
polymer (left). The expansion point is d = a = 4 and the SR exponents hold until the line S. The 
exponent ( is for d fixed a-independent in the SR-regime and a-dependent in the LR-regime, see the 
sketch on the right. 



4.6.1 General structure of the RG flow 

Here, the situation is analyzed as follows. Denote by gi^ and the couplings for SR- 
disorder, and by 6l and 6t the couplings for LR-disorder. Their canonical dimensions are 
in the case of isotropic manifolds 

^l] = [^t] = £ = 2 + L> - , [&l] = [&t] ^S^2 + D- (4.46) 

and for directed manifolds 

k] = bx] = £ = 2 - , = [&t] = 5 = 2- . (4.47) 

Since a < d, we always have S > e. In Sec. 3.10 we have calculated the renormalization 
of elasticity, friction and thermal disorder. To generalize this to the mixed case discussed 
here, one simply has to replace gi^ by gi^ + 6l and g^ by g^ + 6t in the corresponding 
renormalization functions. 

In addition, we now have to take care of the mixture of both types of disorder under 
renormalization. As discussed in Sec. 3.4, no LR-disorder is generated. If we denote 
SR-disorder by — , and LR-disorder by (for isotropic disorder for simplic- 

ity, generalization to transversal and longitudinal disorder are straightforward), then we 
obtain the following egahty for the residues of the pole-terms 




5-1 



Note that this is a statement about the leading term, i.e. the residue, but that the finite 
terms are different and will lead to a less trivial result at the 2-loop level. A similar 
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statement holds for the particle case, as one can see by a procedure analogously to that 
presented in section 3.9. The generalized 1-loop /3-functions then read 

-^B(^L + M(5T + M (4.49) 
Pt{9l, 9t, &t) = -£9t - ~^^9t{9l + &l) + ~ ^) -y-^^T(^T + &t) 

-^'B(^l+&l)(^t + &t) (4.50) 

Pi{9L, 9T, 6l, 6t) = -5 6l - ^Cb^gj^ + 6l) + (l - ^) ^^&l(^t + 6t) (4.51) 

P't{9u 9t: &l, ^t) = -5 6t - ^C6t(^l + &l) + " ^) ^>^^t(^t + &t) (4.52) 

Once fixed points have been found, the exponents can then be obtained, since one has 
simply: 

((fi-L, 9T: ^l, &t) = C{9l + 9t + &t) (4.53) 

and similarly for the exponents z and /3, given by formula (4.8). 

Let us now discuss the structure of these flow-equations. This is not a trivial task 
since the whole depends on the values of the diagrams A, B and C, all of them functions 
of the expansion point, and on S and e. We can furthermore make the observation that 

k:=-^ (4.54) 

is again unchanged under renormalization, thus entering as an additional parameter. The 
situations can then be analyzed by the following observations: 

(1) Relevance of the disorder 

If both e and S are negative, then disorder is indeed irrelevant, i.e. when starting with 
any combination of the couplings, the flow is always towards 0. 

(2) Attractivity of the isotropic fixed point 9 '■= 9-L = 9t 

The best way to study this is to write down the flow equation for the combination 
(fi'L — 9T)/i9L + 9t), which using Eqs. (4.49) and (4.50) cvahiatcs to 



d_ 
^ djj, 



9L-9T {9l - 9t){9l + bL){9T + &t) ^ 
o9L + gT {9l + 9t)^ 



Since the r.h.s. is always positive, the ratio {gi, — 5't)/(5'l + 5't) will flow to 0. We can 
therefore study the simplifled flow-equations for which we set 

9 '■= 9l = 9t 

b := 6t (4.56) 
b]^ — Kb . 
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These equations read 



(5^{g, h) = -eg - ^Cg{g + Kb) + ^Ag{g + 6) - l-B{g + nh){g + h) (4.57) 



2(i •^'•^ ' \ dj 2 "'•^ '2 
P\g, h)^-5h- ^Cb{g + Kb) + (l - ^) + b) (4.58) 

(3) Stability of the SR-fixed point 

Next discuss the stability of the SR-fixed point with 6 = 0, present for e > and 
consequently also 6 > 0. Its stability on the axis 6 = was established before for all cases 
of interest. To analyze its stability towards perturbations of LR-disorder, we linearize Eq. 
(4.58) for small b. Two cases can be distinguished. If (with the diagrams explicitly given 
in Eqs. (4.9) and (4.10)) 

c {d-l){d + 2)A-{d-2)C 

{d-l){d + 2)A-dB-{d-2)C ' ^ ^ ^ 

then LR-disorder is irrelevant and the SR-fixed point alone is stable. The exponents are 
then given by the results for the SR disorder problem. Note that the coefficient on the 
r.h.s. of Eq. (4.59) is always larger than 1, such that this condition is satisfied for weakly 
LR-disorder. 5 < £ would be the result of a naive dimensional estimate, and the above 
results gives the correct result to lowest order, i.e. the tangent to the curve in the a, d 
plane above which SR exponents are correct. For a directed polymer D — 1 one finds the 
condition 

4-a<{4-d)/{l-^) = L1187(4 - d) (4.60) 

OTT 

and for an isotropic polymer 

6 - a < L07723(6 - d) (4.61) 

for the exponents to be given by their SR expressions (see sections 4.2 and 4.4). 

(4) Repulsiveness of the line g = 

By analyzing the fiow-equation Eq. (4.57) for vanishing coupling g, we see that the 
fiow on the axes g = is always pointing upwards (increase of b) at least for «; > 0. 

(5) The complete fiow-diagram 

Let us now draw the complete fiow-diagram for both couplings b and g in the situation 
where the SR-fixed point is stable and k, large. (This domain includes all values of k > 1, 
but is significantly larger. We will detail on the exact phase boundary below.) By direct 
integration of the flow-equations, one finds the situation A4 sketched on Fig. 4.7. The 
figure contains 4 fixed points, namely G, M, C and SR, as expected on general grounds. 
(For directed manifolds, C however is rejected to oo.) How does then the SR-fixed point 
become unstable? 

(6) Crossover from the SR-fixed point to the mixed fixed point M 

Above, we have given the conditions when the SR-fixed point is stable. If one now 
increases 5 such that Eq. (4.59) is no longer satisfied, then the fiow-diagram changes to 



44 





Figure 4.7: Schematic sketch of the possible flow-diagrams for e > 0. For directed manifolds, M and M 
are sufficient and the fixed point C is rejected to infinity. M is relevant when SR is attractive, and iV 
else. The other situations become important for the isotropic case, see the discussion below. 



M on Fig. 4.7. Note that in this case, the mixed fixed point M takes over the stabihty of 
the SR-fixed point. 

(7) Completeness of the scenarios when there is a repulsive SR-fixed point (s > and Eq. 
(4.59) not satisfied). 

We claim that under these conditions, the scenario on the mid of Fig. 4.7 is the only 
possible one, i.e. that no new stable fixed point can appear in the physically relevant 
sector h,g > 0, nor can the fixed point M disappear other than in a reverse of what we 
discussed in point (6), making the SR-fixed point stable. Let us focus on the fixed point 
M. It can neither pass through the axes 6 = other than through the SR-fixed point, as 
else there were 3 fixed points on the axes. Nor can it pass through the axes g = 0, since 
this axes is always repulsive and thus docs not admit a completely stable fixed point. M 
can not pass through the Gaussian fixed point G since 6 > e > 0. The last way for M to 
escape is via infinity. There are situations were this happens (see e.g. [58]), but not here: 
At the fixed point, both Eq. (4.57) and Eq. (4.58) have to vanish. Combining them, one 
obtains for 6 > and g > 

1 (9+^c6)to + ^ 1 

2 g 2 ^ ^ ^ 

bounding g. 
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0.02 0.04 0.06 0.08 ^ 0.1 0.1 0.2 ^ 

Figure 4.8: Some flow-diagrams in the case of isotropic disorder. The left one is for Dq = 1, e = — 1, 
S = 1, K= 100, the right one for Dq = 1, e = 1, S = 2, k = 10. Note also the repulsive SR-fixed point on 
the axis 6 = 0. 

(8) The case of a completely stable SR fixed point 

The same scenario is valid, if £ > and the SR fixed point is attractive (see left of Fig. 
4.7), with two possible exceptions: First, the fixed points C and M can coalesce and then 
annihilate, leaving only two fixed points, G and SR, (flow Q on Fig. 4.7), which always 
have to exist as discussed above. Formally, C and M have become complex. This is only 
possible in the isotropic case, see the discussion below. 

Further changing the parameters, the pair of them can re-appear in the physical sector 
b. g > 0. One then obtains two completely stable fixed points, SR and M, with their 
domains of attraction divided by a phase separatrix, passing through C. The most obvious 
such situation is when k — 0, since then there has to be a fixed point on the axes g — 0, 
see TZ. This is the limiting case of the situation O on Fig. 4.7 for k — > 0. (The axes g — 
is no longer repulsive.) 

(9) The case e < 0, S > 

The last case of interest is, when by say keeping 6 fixed e becomes negative. Then the 
SR-fixed point disappears and only the mixed fixed point stays in the physically relevant 
sector 6 > and g >0. With the same reasoning as above, one concludes that it cannot 
escape, as long as 5 > 0. If however 5 changes sign, it will disappear through the Gaussian 
fixed point at the origin, which then will become stable, as expected. We have plotted 
some of the flow-diagrams in the case e < 0, 5 > on Fig. 4.8, which show amusing spiral 
behavior. We note, as will be discussed below in details that in that e/6 becomes 

very negative, the value of the couplings can become rather large. This terminates the 
discussion of the topology of the fiow-diagrams. 

Note also that in the directed case, the fixed point G is rejected to oo, and thus the 
fiows A4 and Af are sufficient to descrive the crossover. 
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We will now give detailed results for the critical exponents at M and discuss how they 
differ from their values at the LR fixed point discussed earlier. In general, if e and 5 are 
comparable, the critical exponents are continuous functions of both arguments, modifying 
the pure LR-behavior. On the other hand, when 8 is small compared to the modulus of e 
(which then necessarily is negative) then the exponents tend to their values given by the 
LR-fixed points, except in one situation discussed below. 



4.6.2 Results for directed manifolds 

Let us introduce the ratio S = (1 — + 2)A/B. For polymers one expands around 
a = d = 4 and S = 3tt. Then the fixed point M has coordinates: 

bl + g* = Kb*T^ + g* = k{b*^ + g*) (4.64) 

^ = ^TT^hrr (4-65) 
s s 

and is physical, i.e. has a positive = {b^ + g*){l — ^ — + — |) as long as the 
above condition < 1 — ^ is satisfied (it collides with the SR fixed point at this value, 
except for n = 0). The fixed point C in the present case is rejected to 6t = — oo. 

For e/S < 1 — I one thus recovers the previous results for the LR fixed fine of directed 
manifolds, section 4.3, except that k has to be replaced in all expressions (e.g. in Eqs. 
(4.28) - (4.30)) by the value k which depends on the ratio e/S. Thus we have that 
zm{i^) = -2Li?(^) etc.. When e/6 = 1 — ^ one has that all initial k, yield k = 1 (isotropic 
disorder) and one recovers the SR exponents, while when e/6 decreases towards — oo one 
recovers progressively the LR exponents since R, ^ k when e/5 — oo. An interesting 
property is that k{k) is an increasing function which is bounded k{k — oo) — S{1 — 
Thus only a finite interval of the vertical fixed line L studied in the LR expansion can 
really be reached for a fixed finite e/6. 

In particular the Flory exponent for ^ = = (4 — -D)/(2 + a) is still exact at M. This 
yields the exact relation for a* {d) at which the SR exponent will start being corrected as 
a*{d) = (4 — D)/(sji{d) — 2, a relation which could be checked in numerical simulations. 
For polymers using the improved estimates it gives a*{d = 3) ~ 2.76 (while the naive 
extrapolation in Eq. (4.60) gives a*{d = 3) = 2.88) and a*{d — 2) 1.44 (while the naive 
extrapolation in Eq. (4.60) gives a*{d = 3) = 1.76). 

The case k — O must be discussed separately. For concreteness, start with both 6l = 
and gi, — 0. We know that only SR potential disorder is generated, yielding a finite 
g^L, but that 6l remains zero. The short-range part still converges to the isotropic line 
9l — gT = g, as shown in Eq. (4.55) and the above fiow equations are still valid setting 
K — 0. There will be the SR fixed point as before, and the fixed point M is on the axis 
g — 0. At the critical value of e/6 (see Eq. (4.62)) the fixed points SR and M exchange 
stability, but without colliding (right at the critical value there is a fixed line extending 
from one to the other). In any case, the values of the exponents will be as given above 
setting K — 0. 
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Figure 4.9: Domains for which the argument of the root in Eq. (4.68) is negative for polymers (D = 1), 
and therefore the fixed points M and C do not exist. The right figure is just a magnification of the left 
one, on which e.g. domain O cannot bo identified. The dotted line is the boimdary of the domain of 
stability for the SR fixed point as given by Eq. (4.59), such that to the left of it, the SR fixed point 
is stable. Note that it is tangential to the domain, where the fixed points M and C do not exist. We 
have also labeled the regions according to the flow-diagrams in Fig. 4.7. Situation TZ, not drawn here, is 
obtained by starting from the domain O and taking the limit k — » 0. 

We have also checked that the results for the crossover from LR to SR for the particle 
given e.g. in Eq. (62) of [42] are recovered from our formulae. 

4.6.3 Results for isotropic manifolds 

With the same notations as previously (but different coefficients A, B and C) we find the 
coordinates of the fixed points M and C as: 

g* + ly^ = {g* + h*T,Y'' X (4.66) 
g* + hl = g* + Kh*^ = {g* + b*^ f'' K,{X - 1) . (4.67) 

X satisfies an equation with the two roots 

X^^^{l + pG± ^(1 + pGf - 4G) , (4.68) 

where we have defined p— 1 — and G — S{1 — — k). We shall concentrate on the 
case of the polymer D — 1, where one has k,c — 35.8702 and S — ^Kc- We found that in 
the range of interest the fixed point M corresponds to X^ for k > 1 and to X_ for < 1 
and is in the physical sector as long as the condition e/6 < 1 — holds. The fixed 

point C is mostly in the unphysical region. (For the exception sec below.) 

There are several interesting properties. First, one can now study in more detail the 
surprising property of the separatrix found before. One finds that as e/S is decreased 
to — oo at fixed k, two radically different behaviours occur with a transition ai, k — Kc- 
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Either k < Kc and then one finds (for k > 1) that as e/S — oo, saturates at the 
value: 

x+ ^ = (Am) 

while for k > Kc, X_^. keeps on growing indefinitely as | — > — oo as 

In the first case k < Kc the exponents at the fixed point M will evolve smoothly, 
as e/5 ^ — oo to their values on the fixed line obtained in section (4.5) in a first order 
expansion in S at fixed d. However for k > Kc there is no well defined expansion in S at 
fixed d, as the fixed point M becomes non-perturbative in 6 at fixed e < 0. The fixed point 
can still be reached and exponents be computed but only in a perturbative expansion in 
negative e (for 5 > 0). 

For a given k and a fixed ratio | one can compute X using (4.68) and obtain the 
critical exponents at M to first order in S — 2 + D — {2 — D)^ as: 

^=^+(^('^-^)^ ^'-''^ 
z = 2 + ^{R-l){X-l) (4.72) 

/3 = -^(X-l) (4.73) 

with k = Kc(l -1/X), Ko = {d-l)R, R = A/C, and one should set d = 2{2 + D)/{2- D), 
i.e. c? = 6 for the polymer. Using the above limit (4.69) one checks that indeed the 
formulae (4.45) are recovered when e/5 ^ — oo for k < Kc- 

We can now obtain the exponents of the regime k > Kc (which was non-perturbative 
in 5) in an expansion in negative e. We find, using the above limits (4.70): 

C=^ lt£) (IL-l)s (4.74) 

^-^+ (,-'^)7ii) (^-^(^-')^ (^-^5) 

The general analysis of the RG fiow is comphcated and we have therefore restricted 
ourselves to the case of the polymer (D — 1, d ^ dc — 6). Interestingly, there is a domain 
of values of the parameters e/S and k where the two roots given above become complex 
(no physical fixed point M or C). This happens in general for k < 1 when: 

i i — 1^ ^ e 1 I — K 

<l-7<o7^ ^ (4-77) 



5(1+ /^)2 5 5(1- /X)2 

^ y Kc' ^ V Kc' 
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with S = 2/tc/5. For clarity the different regions in the plane of parameters n and 1 — f are 
drawn on Fig. 4.9. The corresponding flow-diagrams are given in Fig. 4.7. We invite the 
reader to start with flow M. and then to continue to AA, C, V, and back to M., and thus to 
observe how the flow is modified. The domain Q is defined by (4.77) and corresponds to 
the case where the two fixed points M and C have annihilated and only the fixed points 
G and SR remain. Note also that the domain O, where l:)oth the SR and the mixed 
fixed points are attractive, and where the both of them are separated by a separatrix as 
shown on O of Fig. 4.7, is very small. Although it is a rather narrow region which may 
thus be difficult to observe numerically (in addition the fiow towards both fixed points 
along the line joining them and passing through C is very slow) it is intriguing physically 
since in that region there will be two different phases and possible values of the exponent 
depending on the respective bare values of LR and SR disorder. 



4.7 Inclusion of self-avoidance 

Up to now we have not taken into account self- avoidance. While this is justified for 
directed polymers (and membranes) by construction, physically, self-avoidance is always 
present for isotropic polymers (and membranes). Although it is difficult to study in a 
controlled way for the present problem, we give here some arguments which indicate that 
self-avoidance docs probably not change the behavior in the new glassy phase. We will 
only address the case of short-range disorder. 

Let us examine the operator for self-avoidance. In the dynamical theory of section 3 
it reads: 

Power counting at the Gaussian fixed point shows that in the absence of disorder, self- 
avoidance is relevant for d < dgj^{D) = AD/ {2 — D). On the other hand we have found 
that short-range disorder alone becomes relevant, for isotropic manifolds, for d < dc{D) — 
(4 -|- 2D)/{2 — D). Since wc study D < 2, self-avoidance is always irrelevant near the 
upper critical dimension dc{D) around which we expand. Thus the calculations of the 
previous sections are consistent. 

There is always the possibility that self-avoidance becomes relevant as the dimension 
d is lowered further. One can give some qualitative estimate. We first consider the naive 
dimension of the self-avoiding operator at the disorder dominated fixed point studied in 
the previous sections. The power-counting used here and in our previous work [48] gives 
for the eigenvalue of the above self-avoiding operator: 

Xl2^D + 2 + (3-{d + 2)C 
where ^F^T-vioiation nieasures the violation of the FDT, 

^FDT-violation ^^_^_^^^_2_/3 + 2C. (4.80) 

Of course, ^F^T-vioiation ^ero in the purely potential case, and one recovers the power- 
counting of self-avoiding membranes without disorder. At the non-trivial fixed point, [3 
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is very small and ( much larger than its Gaussian value 



2-D 



such that (jFDT-violation 



is 



positive and renders self-avoidance less relevant. 

Let us now try to estimate (4.79). If one first inserts the Flory values ^ — 4:/{d + 2) 
and P — 0, one finds that Ag^ — D — 2 + P<0 (since always P < 0) and thus that 
self- avoidance is never relevant. More generally, neglecting /3 (which is always small) one 



finds that self-avoidance becomes relevant if Q\is < 



2+D 
2+d 



The upper bound coincides with 



the Flory value for the roughness exponent in presence of self-avoidance, which is known 
to be a good approximation of the true value. Thus we obtain the rough estimate that 
as long as Cdis > Csa where Csa is the roughness exponent with self- avoidance only, self- 
avoidance can safely be neglected. This is in agreement with the naive expectation, that 
always the operator which results into the larger exponent ( will be dominant. A look 
at the estimates for ^^is obtained in section 4.4 shows that for polymers the neglect of 
self- avoidance should be justified. 

The above approximation is rather naive when applied to the non-trivial fixed point, 
and there are situations where this kind of arguments goes wrong. Since a complete 
analysis of the situation away from £ = is technically impossible, one can only study the 
stability of the fixed point with respect to perturbations by self-avoidance; at first order 
in £ this stability- analysis is then exact^. The complete expression for Asa is 



xfuU 



2D-dC- 5 



FDT-violation 



d 



\nZ. 



p 1 



(4.81) 



flL=9T=9* 



where Zp is the renormalization group Z-factor for self-avoidance. Let us only state the 
result here — • denotes the self-avoidance interaction, see Eq. (4.78)): 



(4.82) 




da — 2 da dr 

In analogy to Eq. (3.78), the diagram is evaluated as 

2D 
2-D 



[C(x,r)+C(y,a)]^-«/^ 



da-2{2- Df 



Since 



d 



InZn 



X ^,,((2 - D)SDr/^ 



>0, 



(4.83) 



(4.84) 



9h=9T=9 

and with the help of Eqs. (4.79) and (4.81), A|^' > A|^, such that vertex-corrections render 
self- avoidance more relevant. The infiuence of these corrections for large e is difficult to 
estimate. Let us only remark, that for large d, {z:xC2^.\^-^ — is exponentially suppressed. 

Thus we obtain that at the isotropic fixed point self- avoidance is presumably not 
important, since the random fiow has already strongly stretched the polymer or membrane. 

^It is worth mentioning that in the case of a 2-body self-avoiding interaction only, 3-body self-avoidance 
is seemingly relevant by power-counting, but that with the help of the above mentioned stability-analysis 
one can falsify this assertion [47]. 
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There is one case however, well studied previously [55,57,59-61], where we know that 
self- avoidance is important. This is when the flow is the gradient of a short-range random 
potential, a case not studied here, but which we mention for completeness. There disorder 
and self- avoidance become relevant simultaneously at d — AD/ (2 — D). The available RG 
treatments for this problem [55,59,60] exhibit a runaway flow to strong coupling, as in 
the present method, and thus the problem can not completely be understood. One open 
question is whether the exponent ( is the same as for pure self-avoidance ^sa [59], and 
even how universal the strong coupling regime is. The question was addressed in [57] 
and it was concluded that the problem is governed by a new strong disorder fixed point, 
very much like the directed polymer problem, where both self-avoidance and disorder are 
relevant. Flory arguments were presented leading to ( ^ 0.8 in c? = 2 in good agreement 
with later simulations [61]. This problem however deserves further studies. 

4.8 Polymer in a hydrodynamic flow: toy model for generation of barriers 

One of the remarkable results obtained in previous sections is that for a manifold with 
D > in a purely divergenceless flow new terms which violate the divergence free condition 
are generated in perturbation theory (this is shown explicitly in Appendix C). Although 
they are flnite, they eventually drive the system to the fixed point where gr = 9l- This 
is in sharp contrast to the particle case D = where the divergenceless flow is a physical 
flxed point, distinct from the flxed point at Qt = q-l- Physically this is because barriers 
are generated in the case of the manifold but not in the case of the particle. 

In the case of the particle the fact that the line g'T = is exactly preserved to all 
orders can be understood by noticing that the exact stationary distribution is known, 
and is simply the spatially uniform distribution. This exact property implies strong 
constraints, as it must exactly be preserved by rcnormalization. In the case of manifolds, 
one can check whether the measure /Xeq = exp(— c/(2T) /^^(Vrj;)^) is also stationary. If it 
was, it would imply similar properties as for the particle, and in particular no generation 
of barriers (since it is translationally invariant). There are some simple cases of flows 
where /icq is indeed stationary. One example is the generalization to polymers [62] of the 
de Marsilly Matheron model of a layered random flow. There one can show explicitly that 
/ieq satisfles the corresponding Fokker Planck equation and is stationary. This problem 
was studied analytically in [62]. Although it can be solved, it is in a sense misleading 
because it misses the physics of the generation of barriers unveiled here: for a polymer 
in a generic divergenceless random flow such a simple translationally invariant stationary 
measure does not exist. 

In this section we will analyze this mechanism further. We will present some toy models 
which allow to illustrate in a simple way how elasticity leads to dynamical generation of 
barriers even in a divergenceless flow. 

We start with a single particle convected in such a flow, which we denote as Vq^ (■^) • 
The probability distribution P(a;, t) for the particle to be at position x at time t, satisfles 
the Fokker Planck equation (T is the temperature): 

dtP{x, t) = d^{Td^P{x, t) - v^{x)P{x, t)) (4.85) 

Since daVa — one checks easily that the spatially uniform measure P{x) — 1/V {V is 
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Figure 4.10: A dumbbell diffusing with each end on a rotating plane, according to Eq. (4.86). 

the volume of the system) is a stationary solution. At large time (and T > 0) the particle 
density will be uniformly distributed over the system. Upon applying a force /, the drift 
velocity in response to the force will be exactly /. This can be checked by methods similar 
to the one introduced in [63]. 

Let us now consider two elastically coupled particles (a dumbbell, see figure 4.10). We 
will start with two particles each located in a 2d plane. For convenience we will denote 
Zk — + Wk the complex position of particle /c = 1, 2 in its plane. We will model the 
directed polymer situation, and thus each particle will see a different flow. The equation 
of motion is: 

ii = c(z2 - zi) + ioJiZi + rji 

Z2 = c{zi — Z2) + ioJ2Z2 + r}2 , (4.86) 

where c is the elastic coupling and = 4T5jj the thermal noise. The flow in each 

plane is thus simply — —uJkV, Vy = uJkX, i.e circulating around the center x — y — 
(where the stream function ip{x, y) — u!k{x^ + y^) has a maximum). The two flows have 
identical center, but a slightly different rotation rate. We denote by 2S = ui — 002 the 
difference, and by ou the mean. Although this model is linear, it already exhibits striking 
behavior. 

In the absence of "disorder", S — 0, the motion is simply a collective rotation around 
the center z — 0. The eigenvalues of the hnear system at T = are iu;, which describes 
the global rotation, and —2c + iuj which describes the decay towards the steady state 
Zi = Z2- At any T the global rotation can be eliminated by redefining Zk = exp{iuit) Zk, 
and similarly redefining the thermal noise. Thus the steady state measure is simply: 

P(zi,Z2) ^e-^\''-'^\\ (4.87) 

i.e. uniform in space and proportional to the Boltzmann measure of the elastic system 
in the rotating frame. Thus again, at T > 0, one observes free collective diffusion of the 
center of mass, and a linear response v = f oi its velocity v to an applied force / . 

Remarkably, as soon as 5 > 0, the zero mode disappears and the two particles converge 
towards the center. At T = the linear system eigenvalues are now A± = — c + io; ± 
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^ (p- — (5^. The global rotation can be eliminated as before, but the differential rotation 
cannot. The convergence of the dumbbell towards [z\ = 0,Z2 = 0) is exponential as 
~ g-(c-Vc2-52)t^ with a characteristic relaxation time r ~ 2c/ at small 5 (which 
interestingly is much larger than the characteristic time 1/\uji—uj2\). The slowest decaying 
eigenmode has Zi/z2 ~ e^^^'^ at small S, i.e. it is twisted by an angle ~ S/c. 

Most interestingly the effect persists to T > 0. It is difficult to find the full stationary 
measure for the process (which seems to be non Gaussian) but it is possible to compute 
explicitly the mean squared positions, by integration of the equations of motion. We find: 

^^l'> = ^(^' + , A; = 1,2. (4.88) 

This indicates a genuine bound state at 2; = 0. At small S, we find l^;^^ ~ Tr with the 
same relaxation time r introduced above. Since T is the diffusion coefficient, the picture 
is that the dumbbell will diffuse away from the center at z = for a time ~ r until it 
is pulled back in again. The nature of this bound state is intriguing. One can easily see 
that the stationary conformation is twisted. Indeed a similar calculation gives: 

{ziz;} = 2T{c + i6)/S\ (4.89) 

Thus, remarkably, as soon as one considers two elastically coupled particles, preferred 
regions appear and the stationary measure is no longer spatially uniform. Physically this 
is because the two particles prefer regions where the elastic energy is reduced. The center 
of mass motion thus becomes coupled to the internal elastic mode. 

Let us now apply an external force. One easily sees that this only amounts to translate 
the center of the bound state. In shifted coordinates the motion remains identical. The 
new center has coordinates: 

a; =F 5 + 2ic , , 

zi,2 = f- 2 x2^ ' 4.90 

where f = fx + ify is the force applied on both monomers. Interestingly the new center 
now depends explicitly on uo (the global rotation cannot be eliminated anymore), in a 
complicated way. Also one has zi ^ Z2i i.e. the dumbbell is slightly inclined. For small 
5 ^ a;, c it reads: 

z,,2^fi—T . ^ . ^ . (4.91) 
\iuj 2cu} — luj^ J 

In the limit 5 — the dumbbell center is close to zi = Z2 = —1/ioJ (with a weak restoring 
spring). Amusingly, in that limit the dumbbell thus responds perpendicularly to the force, 
as a vortex. 

So far we have described situations where the location of the extrema of the stream 
functions in each plane coincide, but the rotation rate does not. If these locations are 
shifted, the same phenomena arise, i.e. a bound state occurs around a new center with 
non trivial position. 

This model can be extended in several directions. First one may consider a directed 
chain with N monomers, again with aligned stream function extrema at 2;^ = 0, and 
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different frequencies One then sees that the decay towards the center at T = exactly 
maps onto the problem of spin depolarization studied in [64], a problem which arises in 
the context of the so called "motional narrowing" in nuclear magnetic resonance (NMR). 
Using the results of [64] leads to a stretched exponential decay of towards the 

center. 

Nonlinear extensions of Eq. (4.86) can also be studied. They will show the generation 
of barriers. For instance in a model with two uncorrelated random stream functions 
ipki^k)-, local long living bound states will be located where extrema are close together 
(note that saddle points are rapidly escaped from) . The dumbbell will eventually escape 
over barriers between these metastable states. 

Finally one can see that the effect persists even if the two monomers see the same 
flow, though in that case it must be nonlinear. A related example of that phenomena 
was discovered by Thual and Fauve [65] in a different context. If one identifies their 
complex Landau Ginsburg order parameter W{x) — z{x) as the complex position in 
2d of monomers in a gaussian elastic chain at T = 0, their complex Landau Ginsburg 
equation describes the dynamics of such a chain in the complex 2d velocity field v = 
{/I + {(3r + i(3i)\z\'^ + {jr + (which describes a double well potential force plus a 

circulating force with a rotation depending on the distance to the center) . They found that 
this equation can have stationary localized solutions, which in our language correspond 
to a rotating and twisted polymer in a bound state. 

To conclude, some insight can be gained by studying these simple models, and explor- 
ing further the physics of these problems would be of great interest. 

5 Conclusion 

In this article, we have studied a model for polymers, and manifolds of internal dimension 
D, diffusing and convected by static random flows. Our results are obtained using fleld 
theoretical dynamical RG near the upper critical dimension dc{D), which is rfc = 4 and 
dc = 6 for directed and isotropic polymers respectively, in short-range correlated flows. 
In principle arbitrary D can be considered, although technically, since dc{D) diverges as 
— > 2~, our quantitative analysis is valid for D < 2. 

We found that below the upper critical dimension dc{D) the random flow is relevant 
and changes the large scale behavior. This yields a new universal RG flxed point for 
polymers and manifolds in non-potential static random flows with short-range correla- 
tions. We have analyzed the critical theory using multilocal operator product expansion 
techniques. As a result, all divergences can be absorbed by multilocal counter-terms, of 
which we explicitly give the 1-loop contribution. 

The new RG flxed points are in general characterized by three exponents, the roughness 
of the manifold C,, the dynamical exponent z, and related to the behavior of the drift 
velocity v{f) under a small applied force /. They have been computed in a dimensional 
expansion, and estimated numerically in physical dimensions. For directed manifolds, a 
relation exists between these three exponents. In the limit D — > we recover previously 
known results for the particle. 

Our results show that in static random nonpotential flows polymers and manifolds 
are in a non-trivial steady state with glassy characteristics. In this state, the polymer or 
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membrane is stretched by the flow with a roughness exponent, which is larger than in the 
absence of the flow (e.g. ( ~ 0.63 in d = 3 and ( a; 0.8 in d = 2 for directed polymers in 
SR random Hows instead of Co = 1/2 in the absence of disorder). The internal dynamics is 
slower, with z > 2. Upon applying a small uniform external force /, the polymer will move 
with a reduced velocity v{f) ~ with (p > 1. This state is characterized by a zero linear 
mobility ji = ^|y=o- Thus one can consider that the polymer is trapped in preferential 
regions of the flow. On a simpler toy model with only two particles, we have been able 
to describe more precisely the effect of trapping, but in the case of the polymer and 
manifolds its detailed physics remains to be investigated. Numerical simulations would 
be extremely useful to check our results for the exponents and to further characterize 
this glassy phase. These simulations should be easier to perform than for conventional 
potential glassy systems, since in the present case the convergence to the steady state is 
expected to be faster. 

The glassy state that we found for this non-potential system can be compared to the 
one arising in potential systems. The main difference is that as soon as disorder is non- 
potential, an effective temperature is generated, resulting in a violation of the fluctuation 
dissipation theorem. The increase of temperature competes with the relevance of disorder 
which tends to make the temperature irrelevant. Since the important dimensionless dis- 
order parameter is proportional to the inverse temperature, contrarily to most potential 
systems, disorder docs not flow to strong coupling. Instead, it flows to a perturbatively 
accessible fixed point, studied here. The effect of disorder, and thus the glassy character- 
istics, are in effect weaker than in potential glasses. Indeed, we find a finite dynamical 
exponent z > 2 and f ~ /"^ rather than the much stronger non linear "creep" behavior 
V ~ ex.p{—c{fc/fY) obtained for manifolds with D < A and z effectively infinite, which 
arises from a zero temperature fixed point [12,15,66]. In fact the state found here is remi- 
niscent of the so called "marginal glasses" predicted in the statics for periodic manifolds at 
their lower critical dimension D = 2 [49,67-72], where also thermal effects are important 
and which have a similar v ^ characteristics. This corresponds effectively to barriers 
which grow logarithmically with size (and inverse applied force), and it is remarkable 
that the "randomly driven polymer" also exhibits such divergent barriers (in an effective 
way, since the notion of barriers is not clear-cut in nonpotential systems). It is all the 
more remarkable in the case of a purely hydrodynamic (i.e. divcrgcncclcss) flow, which 
for the single particle leads to a trivial, spatially uniform, stationary distribution (thus 
without barriers). Here, these barriers are generated from a non trivial interplay between 
elasticity, disorder and thermal fluctuations. 

In the case of LR correlated flows we have found via an expansion in the range of 
the correlation a near a = ac{D) = dc{D) at fixed d that the behavior of the system is 
controlled by a fine of fixed points, and thus depends continuously on the ratio n = gh/dT 
of the longitudinal to transversal part of disorder (which remains unchanged under renor- 
malization). For directed manifolds, our analytical results were found to be in agreement 
with the exact large d result obtained in [44] and even, to first order in e with the cor- 
responding Hartree approximation. In that case ( was found to be exactly equal to its 
Flory approximation Cf = (4 — -D)/(2 -|- a) while the other exponents are continuously 
varying. Several remarkable results were obtained in the case of the isotropic manifold. 
We found that there exists a critical value Kcr of k, beyond which no fixed point is found 
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in the fixed dimension expansion, as well as a smaller characteristic value Kq beyond 
which the manifold becomes compressed C < Co rather than stretched by the flow (as it 
is in all other cases). Since this problem also corresponds to directed polymers subjected 
to non-potential "correlated disorder" , these effects may be understood as a result of a 
competition between localization (due to the correlated nature of the force) and driving. 
To further investigate this competition we also studied in detail the intricate crossover 
between SR and LR disorder, in a dimensional expansion around the point a = d = 6 
for polymers, important for studies in the physical dimensions d — 2,3. In view of the 
rich variety of obtained behaviors, numerical simulations would be welcome to study the 
competition between SR and LR disorder effects. 

We also note that the glassy phase obtained here is via the RG perturbatively con- 
trolled around the line dc{D) (see Fig. 2.1), whereas the large d approach of [44] works for 
any D and for LR disorder. A last and complementary method is the functional renor- 
malization group (FRG) which can be used in a D — 4 — £ expansion and thus allows to 
explore the region of larger D than the present method. Recent progresses on the finite 
temperature FRG in [33] which showed the existence of such a finite temperature fixed 
point, leads us to conjecture that a similar fixed point will be accessible by FRG in the 
present model and should merge smoothly into the one found here upon decreasing D. 
Work is in progress to verify this conjecture [73] . 

To conclude, let us mention an interesting speculation, inviting for further investigation 
[74]: It is well-known that a directed polymer in a random potential maps via the Cole- 
Hopf transformation to the Burgers-Kardar-Parisi- Zhang equation [16]. The question is 
whether the generalization of the directed polymer diffusing in a non-potential random 
force field instead of a potential one, has some effective correspondence in a generalization 
of the Burgers equation to non-potential driving and is thus possibly related to the Navier 
Stokes equation. 
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A Correlator and Response-function 

In this appendix we give a collection of formulas for the correlator and the response 
function in the Ito-discretization. 

1 9 A 

- {r{k, u)r{k\ J)), = -r^^^, {27rf S^{k + k'){2^)5{u + J) (A.l) 

- {r{k,uj)f{k',u;')), = {2n)^d^ik + k')i2n)d{u; + J) (A.2) 
d luj -\- Xk'^ 

-lf{k,uy{k',uj')), = Q. (A.3) 
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Fourier-transformation with respect to uj yields 
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d 



'Ak,t)r{k',t% 



1 



-A|t-t'|fe2 



d 

^-{f{k,t)f{k',t% = 0. 



(A.4) 
(A.5) 
(A.6) 



We will mostly use the real space representation. Let us start with the correlator. 

The Fourier-transformation of Eq. (A.4) cannot be done explicitly but involves an error- 
function. In addition, to obtain a well-defined function for D < 2 one has to subtract the 
zero-mode. We therefore define 



C{x,t) :=^((r(x,t)-r(0,0)f;^ 



The hmiting cases are 

C{x,t) = 

with 

Let us also note the ratio 
The response-function 
is calculated analytically 
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R{x,t) :=-(r(a;,t)f(0,0))o 



Rix,t) = eit > 0) (47rAi)-''/'e-^'/^^* . 
The (perturbative) fluctuation dissipation theorem (FDT) reads 

e{t)-^Cix,t) = R{x,t) . 



(A.7) 



(A.8) 



(A.9) 
(A.IO) 
(A.11) 

(A.12) 
(A.13) 



B Extrapolations 

The numerical values for the exponents reported in table 4.2 and 4.4 were obtained by 
an optimization procedure using the freedom to expand about any point on the critical 
curve e{D,d) = 0. This procedure, first proposed in [75], was later refined in [52] to 
obtain reliable estimates for the scaling dimension of self-avoiding polymerized tethered 
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membranes in the frame work of a 2-loop calculation. As will become apparent, various 
extrapolation schemes are possible, and choosing the best one is almost an art; we shall 
rely heavily on the methods developed in Ref. [52] to which the interested reader is 
referred to for further details and discussion. Our presentation is inspired by [58], where 
complementary material can be found. The general idea is of course to expand about some 
point {Do, do) on the critical curve e{Do, do) = 0. The simplest scheme is to extrapolate 
towards the physical theories for D = 1,2 and d = 2,3, . . ., using the expansion parameters 
D — Do and d — do. However, as shown in Ref. [52], this set of expansion parameters is 
not optimal, and better results are obtained by using Dc{d), (solution of the equation 
e{Dc{d),d) = 0) and e{D,d). Furthermore, it is advantageous to make expansions for 
quantities such as (d rather than (. 

After selecting one of these schemes, the next step is to re-express the quantity which 
shall be expanded in terms of Dc{d) and e{D,d). If we are interested in polymers {D = 1) 
in d — ?> subject to isotropic disorder, we have to evaluate this expression for e — 3/2. 
However, we are still free to choose the expansion point along the critical curve, which 
then fixes Do- As the expansion point is varied, different values for Cd (i.e. 3^ for d = ?>) 
are obtained, as plotted on Fig. B.l. The criterion for selecting a value for C, from such 
curves is that of minimal sensitivity to the expansion point Do- We thus evaluate C, at 
the extrema of the curves. The broadness of the extremum provides a measure of the 
goodness of the result, and the expansion scheme. The robustness of this choice was 
explicitly checked in Ref. [52] in the case of self-avoiding tethered membranes, by going to 
the second order. (For additional discussions of such "plateau phenomena" see Sec. 12.3 
of Ref. [52].) 

While we examined several such curves, only a selection is reproduced in Figs. B.l, 
B.2, and B.3. While quantitative predictions were already reproduced in tables 4.2 and 
4.4, let us focus on some important qualitative observations. First, it is striking to observe 
that the dependence of on D in the directed case and the independence on D in the 
isotropic case, read off from the extrapolations on Fig. B.2, is already predicted by the 
Flory estimate. Indeed, the latter is hard to improve numerically. 
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Figure B.l: Extrapolations for ^rf to (D = 1, d = 3). 
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Figure B.2: Extrapolations for ^ to polymers and membranes in 3 dimensions, using the extrapolation of 
(d. The left diagram is for the directed, the right for the isotropic case. 



For the exponents z and (3, we did not find a reliable extrapolation scheme: As can be 
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Figure B.3: Extrapolations for /3 (left) and z (right) to polymers and membranes in 3 dimensions, isotropic 
disorder. Note that extrapolations are impossible for z and unreliable for /3. 



deduced from Fig. B.3, the 'minimum sensitivity point' in extrapolations for f3 is roughly 
independent of Dq, which indicates that the extrapolation for /? only depends on e, but 
not in a non-trivial way on D and d. 

Extrapolations for z show a plateau when approaching Dq ^ 2 only. Since the struc- 
ture of the perturbation expansion is such that expanding about Dq — 2 always yields 
the free-theory result = 2, extrapolations merely permit in concluding on z > 2. On 
the other hand, it is possible to build up combinations of z and D and d, which upon 
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expansion build up a nice plateau. However the answer obtained by this method strongly 
depends on the selected scheme, and therefore is unreliable. 



Finite corrections to longitudinal disorder in the purely 

transversal case 



The aim of this appendix is to show that in the pure transversal case, longitudinal disorder 
is generated at the 2-loop order. A full 2-loop calculation is very difficult to perform 
(see [52]). We will therefore only show, that at least one MOPEl-coefficicnt, namely 

-"-^^^^^ does not vanish, as this was the case for t -"^-^^^^ . Since already 

the evaluation of t -^A^^ is a formidable task, we simplify the task further by 

considering only a special limit. This is accomplished by evaluating the leading singularity 
when approaching two of its vertices as 
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The symbol — hereby denotes the subleading operator (with the possible transversal 
or longitudinal projectors unspecified) 



= / f(a;,t)A;2e''=[''(^'*)-''(f'*')]f(2/,0 . 



(C.2) 



Since we know from the first-order calculations, see section 3.4, that 



0, the first two terms on the r.h.s. in Eq. (C.l) vanish. Prom the cal- 

= 0, but that r:^ 



the next-to-leading term in the 



is non-zero m 



culations in section 3.4 also follows that (z 
does not vanish. Since moreover — is after 
MOPE, at least for D — > 2, it is sufficient to show that 

order to establish that ^ ^ does not vanish identically. This is the aim of 

the following calculation. Let us however note that the argument is not conclusive, when 
--^"^ is not the most relevant subleading calculation. This appfies especially to the case 
of the particle with D — 0, where the canonical dimension of Vr becomes and there is 
thus an infinity of marginal operators, which can be constructed from — — ^ by adding 
an arbitray power of Vr. We comment on this special case later, see appendix D. 
We now study the MOPE 



+ derivatives of w>»- 
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To simplify notations, we give all results for the isotropic case and short-range disorder 
{a = d). We start from the expression (3.34) for the contraction Di in the notation of 
figure 3.1 and the analogue expression for D2, which together read 



2 / / ^ik[rix,0)-r(y,t}] g(-pH^)[C(x-.',r)+C(,-,',a)]^(^ _ ^> ^ ^^^^y _ y, ^ ^) ^ 

Jk Jp 

[f«(x, o)f^(i/, t) {p - (p - 1)5 - f«(x, 0) h{y, t) {p - f )^ {p + ^)(3 
[grPlpiP - I) + g^P^^ip - 1)] [gTP^^sip + 1) + g^P^sip + 1 

We want to write this expression in the form 

/ f„(x,0)f^(y,i)e*'=['^(-'°)-'-(^'*)]A,^(A;) . 

J k 



(C.4) 



(C.5) 



In order to decompose Aap{k) into longitudinal and transversal part, we remark that if 
we denote 

A^f^ik) = A^{k)P^^{k) + A^{k)P^^{k) , (C.6) 

then 



AL(fc) = pi:p{k)A^m 

At(A;) = ^PJ^(^)A„^(A;) . 



(C.7) 



Applying the longitudinal projector P^p{k) — onto Aaf3{k) as obtained from Eqs. 
(C.4) and (C.5) yields 
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Using the FDT (3.12) and expanding up to second order in k gives for the longitudinal 
part 

d_d_^^p2[C{x-x',T)+C{y-y',a)] 

p da dr 
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Two remarks are in order: First, integrating over p, we recover from the first term pro- 
portional to QlOt the leading contribution already given in Eq. (3.37). Second, and most 
importantly, there is a term proportional to g^, which is positive for all dimension d > 1. 
Thus, even when starting with pure transversal disorder, longitudinal disorder will be 
generated under renormalization. Although finite at 1-loop order, it will necessitate a 
counter-term at the 2-loop level, as discussed at the beginning of this section. 

For completeness, let us also analyze corrections to the transversal part of disorder. 
Since it is calculatory more convenient, instead of contracting Aq^(/c) with the transversal 
projector, we contract it with the unit-matrix 5^^, resulting in {d — 1) At + Al. It is then 
easy to subtract Al already given in Eqs. (C.8) and (C.9). Ergo 
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Using the FDT, Eq. (3.12), and expanding up to order /c^, this can be written as 
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The transversal part can now be obtained by subtracting Eq. (C.9) from Eq. (C.ll). The 
interesting thing to note is that there is no component proportional to g^ as expected from 
the fluctuation dissipation theorem, which states that potential disorder is preserved under 
renormalization. On the other hand, there will be all other combinations of terms. 



D Non-renormalization of transversal disorder in the particle case 

In this appendix, we show, why in the particle case purely transversal disorder is not 
renormalized. An analogous result in the context of the Focker-Planck equation was 
obtained in [42]. 

To this aim consider n interaction vertices — ^ , and choose from its 2n end-points 
the one who is the most retarded in time. Without loss of generality due to translation- 
invariance in time its time-argument can be set to 0, leading to ra{0)- (Note that for a 
particle, there is only a time argument, but no space argument x.) To construct a new 
interaction with exactly two response-fields per vertex from the n interaction vertices, it 
is necessary to contract fa(0). Now, any of the other vertices has the structure V := 
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j^^(^^)ga[r(t)-r(t')]^(i'). Contracting f„(0) with V yields ik^R{t) - iKR{t'). Since by 
assumption t > and t' > 0, R{t) = R{t') = 1 and the both terms canceL This cancelation 
is independent of the kind of disorder, and also present for longitudinal one. (Note, that 
the argument is not valid for D > 0.) The difference between the two kinds of disorder 
comes from the last possible contraction, namely by contracting fQ,(0) with the other 

rp /T 

end of the same vertex. Denoting by P^^ (/c) the transversal respectively longitudinal 

projector, this vertex is J/. Pj^^(A;)fQ(0)e*'^[''(°)~''Wlf^(t)^ and contracting ra{0) with e"**^''^*) 

yields —ika which is multiplied by P^p^{k). In the transversal case, the product kaPjf^ik) 
vanishes identically. This is not so for longitudinal disorder, leading to a renormalization 
of disorder. 
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